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ments for  the  degree  of  Master  of  Science. 


ABSTRACT 


This  thesis  investigates  the  use  of  cwo  types  of 
stability  analysis,  i.e.  total  stress  and  effective  stress, 
for  evaluating  the  factor  of  safety  of  embankments  constructed 
on  soft  clay.  Three  classes  of  stability  problems  are 
defined  based  upon  the  drainage  conditions  during  construction 
and  shear.  The  use  of  both  types  of  analyses  for  each  class 
of  stability  is  discussed. 

The  factors  of  safety  computed  for  both  types  of  analyses 
are  compared  for  three  situations:  eight  case  studies  of 
actual  embankment  failures;  variation  with  embankment  height 
for  seven  embankments;  and  the  long  term  stability  for  two 
hypothetical  embankments.  ^ 

Two  cases  received  special  attention:  variation  in  the 
factor  of  safety  with  height  for  a test  embankment  constructed 
to  failure  in  Portsmouth,  N.H.  and  the  long  term  stability  of 
a hypothetical  embankment  resting  on  a soil  having  the 
anisotropic  strength  properties  of  Boston  Blue  Clay. 

For  design  of  embankments  and  determination  of  the  end 
of  construction  stability  the  total  stress  analysis  is 
recommended  as  being  more  reliable  and  easier  to  apply.  The 
undrained  strength  should  be  determined  by  applying  Bjerrum's 
empirical  correction  factor  to  field  vane  tests  or  from 
normalized  soil  properties  such  as  the  SHANSEP  technique. 

For  variation  of  the  factor  of  safety  with  embankment 
height,  either  the  total  or  effective  stress  analysis  may  be 
used  but  the  effective  stress  analysis  will  always  give  a 
higher  factor  of  safety  for  the  same  height  of  embankment. 

For  long  term  stability  the  effective  stress  analysis 
is  easy  to  apply  and  will  give  results  significantly  higher 
than  a total  stress  analysis  but  in  either  case  the  long  term 
situation  is  not  the  critical  class  of  embankment  stability. 


Thesis  Supervisor: 
Title : 


Charles  C.  Ladd 
Professor  of  Civil  Engineering 
-2- 


r\ 


ACKNOWLEDGMENTS 


I would  like  to  acknowledge  Professor  Charles  C.  Ladd, 
my  thesis  supervisor,  for  his  teaching,  help  and  encouragement 
throughout  my  stay  at  M.I.T. 

The  United  States  Air  Force  is  recognized  for  their 
interest  and  financial  support  of  my  graduate  education. 

A special  thanks  to  my  wife  who  gave  me  the  encouragement 
to  complete  my  studies  and  to  my  daughter  who  deserved  the 
time. 


NTIS  GRAJfcl 


DDC  TAB 

Unannounced 

Justification 


ribution 


Awi 7 a M. llt.y  Codes 
| Avail  and/or 


TABLE  OF  CONTENTS 


i 

1 


. I 


TITLE  PAGE 

ABSTRACT 

ACKNOWLEDGMENTS 

TABLE  OF  CONTENTS 

LIST  OF  TABLES 

LIST  OF  FIGURES 

LIST  OF  SYMBOLS 

CHAPTER  1 INTRODUCTION 

CHAPTER  2 METHODS  OF  STABILITY  ANALYSIS 

CHAPTER  3 TYPES  OF  ANALYSES 

3.1  TOTAL  STRESS  ANALYSIS 

3.1- 1  UNDRAINED  SOIL  BEHAVIOR 

3.1- 2  DETERMINATION  OF  UNDRAINED  STRENGTH 

3.2  EFFECTIVE  STRESS  ANALYSIS 

CHAPTER  4 CLASSES  OF  STABILITY 

4.1  "UU"  END-OF-CONSTRUCTION 

4.2  "CD"  LONG  TERM 

4.3  "CU"  INTERMEDIATE 

CHAPTER  5 COMPARISON  OF  TOTAL  AND  EFFECTIVE  STRESS 

ANALYSES 

5.1  CASES  OF  EMBANKMENT  FAILURE 

5.2  VARIATION  WITH  EMBANKMENT  HEIGHT 

-4- 


-«t.  t 


; ■*' 


r - >•  ..  - ' 


PAGE 


1 

2 

3 

4 
6 
7 
9 

13 

17 

22 

23 

24 
26 
29 

35 

36 
40 
42 
49 


50 

52 


TABLE  OF  CONTENTS  CONTINUED 


LONG  TERM  STABILITY 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  SUMMARY  AND  CONCLUSIONS 

6.2  RECOMMENDATIONS 


PAGE 

55 

67 

67 

71 


REFERENCES  73 
APPENDIX  A PORTSMOUTH,  NEW  HAMPSHIRE  TEST  EMBANKMENT  78 
APPENDIX  B MODEL  EMBANKMENT  ON  BOSTON  BLUE  CLAY  92 


LIST  OF  TABLES 


TABLE 


TITLE 


PAGE 


EQUILIBRIUM  CONDITIONS  SATISFIED  BY  VARIOUS 
METHODS  OF  STABILITY  ANALYSIS 

TSA  vs  ESA  FACTOR  OF  SAFETY,  ACTUAL  EMBANKMENT 
FAILURES 


TSA  vs  ESA  FACTOR  OF  SAFETY  VARIATION  WITH 
EMBANKMENT  HEIGHT 

TSA  vs  ESA  FACTOR  OF  SAFETY  MODEL  EMBANKMENTS 


INPUT  GEOMETRY,  PORTSMOUTH  N.H.  EMBANKMENT 


SOIL  PARAMETERS,  PORTSMOUTH  N.H.  EMBANKMENT 

PIEZOMETER  DATA,  PORTSMOUTH  N.H.  EMBANKMENT 

TOTAL  STRESS  ANALYSIS  vs  EMBANKMENT  HEIGHT, 
PORTSMOUTH  N.H. 


EFFECTIVE  STRESS  ANALYSIS  vs  EMBANKMENT  HEIGHT, 
PORTSMOUTH  N.H. 

MANUAL  EFFECTIVE  STRESS  ANALYSIS,  PORTSMOUTH 
N.H.  (I) 

MANUAL  EFFECTIVE  STRESS  ANALYSIS,  PORTSMOUTH 
N.H.  (II) 

UNDRAINED  STRENGTH  RATIOS  CORRECTED  FOR  STRAIN 
COMPATIBILITY  vs  OCR  FOR  BOSTON  BLUE  CLAY 

SETTLEMENTS  AT  THE  CENTERLINE,  BBC  MODEL 
EMBANKMENT 


l 


UNDRAINED  STRENGTH  BEFORE  CONSOLIDATION,  BBC 
MODEL  EMBANKMENT 

UNDRAINED  STRENGTH  AFTER  CONSOLIDATION  FOR 
SLOPE  OF  BBC  MODEL  EMBANKMENT 

UNDRAINED  STRENGTH  AFTER  CONSOLIDATION  FOR 
CREST  OF  BBC  MODEL  EMBANKMENT 


■ ••  • *,  r,  ' * 

w ...o  v .• 


LIST  OF  FIGURES 


.1 


B 


TT 

L 


n 


I 

I 


FIGURE 


TITLE 


PAGE 


3-1 

3- 2 

4- 1 

4-2 

4- 3 

5- 1 


5-2 

5-3 

5-4 


5-5 

A-l 

A-2 


A-3 


B-l 

B-2 

B-3 

B-4 


B-5 

B-6 


« 0 ANALYSIS 
"TYPICAL"  STRESS  PATH 

YLIGHT  MODEL  OF  PORE  PRESSURE  PREDICTION 
STRESS  PATH  FOR  LONG  TERM  STABILITY  ANALYSIS 
FACTOR  OF  SAFETY  vs  EMBANKMENT  HEIGHT 


ACTUAL  EMBANKMENT  FAILURES  TSA  vs  ESA  FACTOR 
OF  SAFETY 


33 

34 

46 

47 

48 

62 


ACTUAL  EMBANKMENT  FAILURES  TSA  vs  ESA  FACTOR 
OF  SAFETY  CORRECTED  FIELD  VANE  AND  SHANSEP 


FACTOR  OF  SAFETY  vs  EMBANKMENT  HEIGHT, 
PORTSMOUTH  N.H. 


TSA  vs  ESA  FACTOR  OF  SAFETY  WITH  VARIATION  IN 
EMBANKMENT  HEIGHT 


63 

64 

65 


CONNECTICUT  VALLEY  VARVED  CLAY  MODEL  EMBANKMENT  66 


SOIL  PROPERTIES,  PORTSMOUTH  N.H. 


89 


CROSS  SECTION  OF  EMBANKMENT  FAILURE,  PORTSMOUTH  90 

N.H. 


STRENGTH  ALONG  FAILURE  SURFACE,  TSA  vs  ESA 
PORTSMOUTH  N.H. 


91 


STRESS  HISTORY  BBC  MODEL  EMBANKMENT 
EMBANKMENT  GEOMETRY  BBC  MODEL  EMBANKMENT 
UNDRAINED  STRENGTH  BOSTON  BLUE  CLAY 


ANISOTROPIC  UNDRAINED  STRENGTH  PARAMETERS 
BOSTON  BLUE  CLAY 


102 

103 

104 

105 


UNDRAINED  STRENGTH  BBC  MODEL  EMBANKMENT 


BEFORE  CONSOLIDATION  FAILURE  SURFACE  AND 
STRENGTH  BBC  MODEL  EMBANKMENT 


106 

10? 


-7- 


' ' • •*;  - >*  ♦ \ A.  ' 


LIST  OF  FIGURES  CONTINUED 


T 


FIGURE  TITLE 

B-7  AFTER  CONSOLIDATION  FAILURE  SURFACE  AND 
STRENGTH  BBC  MODEL  EMBANKMENT 


PAGE 

108 


LIST  OF  SYMBOLS 


SYMBOL 


CK0U 


cvvc 


R 


i(su(v)  + su(h) ) 

Active,  Direct  and  Passive 
(rff(45°))/(Tff(H)  xTff(V} 

Boston  Blue  Clay 
Undrained  strength 

Cohesion  intercept  for  effective  stress 

Anisotropically  Consolidated  Undrained  with  pore 
pressure  measurements 

Consolidated  Drained 

Isotropically  Consolidated  Undrained  with  pore 
pressure  measurements 

Consolidated  K0  Undrained  with  pore  pressure 
measurements 

Compression  Ratio 

Consolidated  Undrained 

Consolidated  Undrained  with  pore  pressure 
measurements 

Coefficient  of  consolidation 
Connecticut  Valley  Varved  Clay 
1 - (b/a)2 
Direct  Simple  Shear 
(1-KS)/(1+KS) 

Effective  Stress  Analysis 
Effective  Stress  Path 


. ' 1 - - # t, ' + » . • 

'r*  V;  • *«  * • 

*.«►  ‘V  - * j/l 


LIST  OF  SYMBOLS  CONTINUED 


SYMBOL 

FS 

FV 

H 

ICES 

Ko 

Ks 

KSF 

L 

LEASE  II 
M 

OCR 

P 

*c 

PCF 

PSC 

PSE 

PSF 

q 

R 

RR 

Sd 

Su 


Factor  of  Safety 
Field  Vane 
Horizontal 

Integrated  Civil  Engineering  System 
Lateral  stress  ratio  cr  (H) / (V ) 

Anisotropic  strength  ratio  (H)/Pff (V) 

Kips  per  Square  Foot 

Loading 

Limiting  Equilibrium  And  Soil  Engineering  II 
Moment 

OverConsolidation  Ratio  ^vm/cvo 
Average  deviator  stress  (<r ^+^3)  /2 
Preconsolidation  load 
Pounds  per  Cubic  Foot 
Plane  Strain  Compression 
Plane  Strain  Extension 
Pounds  per  Square  Foot 
Maximum  shear  stress  /2 

Maximum  shear  stress  at  failure 
Radius  of  failure  circle 
Recompression  Ratio 
Drained  shear  strength 
Undrained  shear  strength 


-10- 


LIST  OF  SYMBOLS  CONTINUED 


SYMBOL 

SHANSEP 

STAB3D 

TSA 

TSP 

u 

u 

e 

uo 

U 

uc 

uu 

V 

Y LIGHT 


Stress  History  And  Normalized  Soil  Engineering 
Properties 

Program  for  3D  stability  analysis,  includes 
anisotropic  strength  capability 

Total  Stress  Analysis 

Total  Stress  Path 

Pore  pressure 

Excess  pore  pressure 

Initial  pore  pressure 

Unloading 

Unconfined  Compression 

Unconsolidated-Undrained 

Vertical 

Yield  Locus  Influenced  by  Geological  History  and 
Time 


*t 

e 

£ cf 

<r 

c 

<t3 


Unit  weight 
Total  unit  weight 

Inclination  of  the  failure  surface 
Bjerrum’s  field  vane  correction  factor 
Final  consolidation  settlement 
Total  normal  stress 
Effective  normal  stress 
Major  principal  stress 
Minor  principal  stress 


-11- 


LIST  OF  SYMBOLS  CONTINUED 


SYMBOL 

Srf 

*v 

^ vc 

3vf 

<J"vm 
3r  vo 

r 

^ ff 

4 

1 


Effective  normal  stress  at  failure 
Effective  vertical  stress 
Effective  vertical  consolidation  stress 
Effective  vertical  final  stress 
Maximum  vertical  effective  stress 
Initial  vertical  effective  stress 
Shear  stress 

Shear  stress  on  the  failure  surface  at  failure 
Friction  angle 

Friction  angle  for  effective  stress 


-12- 


CHAPTER  1 


INTRODUCTION 

This  thesis  discusses  and  compares  the  total  stress 
analysis  (TSA)  and  the  effective  stress  analysis  (ESA)  for 
estimating  the  stability  of  embankments  constructed  on  soft 
clay  foundations.  Although  both  types  of  analysis  have 
their  theoretical  justification,  the  factors  of  safety 
calculated  from  each  produce  different  numerical  results 
and  their  use  and  interpretation  have  caused  much  debate 
amoung  geotechnical  engineers.  This  research  is  an  effort 
to  sort  out  the  advantages  and  disadvantages  of  each  type 
of  analysis  and  make  recommendations  for  their  use. 

"If  we  are  wise  we  will  learn  all  we  can  about  the 
implications  of  both  methods  (TSA  and  ESA)  and  be 
in  a position  to  apply  whichever  one  best  suits  the 
situation  at  hand."  (Peck  and  Lowe,  i960) 

The  study  is  limited  to  embankments  constructed  on 
soft  clays  and  does  not  include  excavations  or  natural 
slopes  in  clays  or  construction  on  granular  materials.  The 
term  "soft  clay"  implies  that  embankment  construction  will 
load  the  clay  foundation  beyond  its  preconsolidation  stress 
and  into  the  normally  consolidated  range.  "Stiff  clays", 
i.e.  those  which  remain  overconsolidated  after  loading 
would  usually  not  have  a problem  with  stability  and 


therefore  are  not  considered.  Despite  this  limitation  in 
scope  the  profession  frequently  encounters  situations 
involving  embankments  constructed  on  soft  clay,  for  example 
in  transportation  embankments  (highways,  railroads,  rapid 
transit  and  bridge  approaches),  levees,  earth  dams  and  for 
material  stockpiles.  Because  embankments  often  involve 
large  volumes  of  earthwork,  it  may  be  economical  to  design 
embankments  with  factors  of  safety  of  1.3  or  lower.  When 
using  such  relatively  low  factors  of  safety,  the  accuracy 
of  the  stability  analysis  becomes  increasingly  important 
for  proper  design. 

Chapters  2 through  4 will  discuss  'methods'  of  stability 
analysis,  the  'types'  of  stability  analysis  and  the 
'classes'  of  stability  repectively.  The  'method'  of 
stability  analysis  involves  the  choice  of  the  potential 
failure  surface  and  the  assumptions  necessary  to  make 
statically  determinate  the  equations  of  equilibrium.  All 
methods  of  analysis  use  a limiting  equilibrium  definition 
for  the  factor  of  safety  which  is  equal  to  the  ratio  of 
the  shear  strength  of  the  soil  along  the  assumed  failure 
surface  to  the  shear  stress  required  for  equilibrium. 
Possible  sources  of  error  in  this  calculation  are  briefly 
discussed. 

The  'type'  of  stability  analysis  refers  to  the  form  of 
the  calculation  used  in  determining  the  shear  strength  of 
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the  soil.  This  expression  of  the  shear  strength  is  based 
either  on  a total  stress  analysis  or  an  effective  stress 
analysis.  In  a total  stress  analysis  the  shear  strength 
along  the  failure  surface  is  assumed  to  be  independent  of 
the  applied  total  stress,  but  of  course  may  vary  with  depth 
(changes  in  stress  history)  and  inclination  of  the  failure 
surface  (due  to  anisotropy).  The  effective  stress  analysis 
assumes  an  expression  for  the  shear  strength  along  the 
failure  surface  which  is  dependent  on  the  effective  normal 
stress,  i.e.  total  stress  minus  pore  pressure.  Procedures 
for  measuring  the  strength  parameters  for  both  types  of 
analyses  are  also  discussed. 

Three  'classes’  of  stability  problems  are  defined 
based  upon  the  drainage  conditions  assumed  with  respect  to 
consolidation  and  shear.  These  classes  are  labeled  after 
similar  laboratory  test  conditions:  "UU",  unconsolidated 
undrained  or  end-of-construction;  "CU" , consolidated 
undrained  or  intermediate;  and  "CD",  consolidated  drained 
or  long  term  stability  analysis.  The  application  of  both 
types  of  stability  analysis  is  discussed  for  each  class  of 
stability. 

Chapter  5 contains  case  studies  and  examples  comparing 
total  stress  and  effective  stress  stability  analyses  of 
embankments  on  soft  clay.  The  first  section  compares  the 
results  of  several  embankment  failures  reported  in  literature 
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which  were  analyzed  by  both  types  of  analysis.  The  second 
section  compares  the  variation  in  the  factor  of  safety  with 
embankment  height  when  computed  by  both  types  of  stability 
analyses  and  the  third  section  compares  the  factor  of  safety 
of  two  fully  consolidated  embankments. 

The  final  chapter  contains  a summary  and  presents 
recommendations  for  the  use  of  total  stress  and  effective 
stress  analyses  for  embankments  on  soft  clay. 


CHAPTER  2 


METHODS  OF  STABILITY  ANALYSIS 

In  any  earth  structure  or  foundation  which  applies 
load  to  the  soil,  a most  important  consideration  for  design 
is  deformations.  To  keep  deformations  within  acceptable 
bounds,  two  separate  calculations  are  usually  performed: 
a stability  analysis  and  a settlement  analysis.  The 
stability  analysis  is  performed  in  an  effort  to  assess  the 
possibility  of  a gross  shear  failure  beneath  the  structure 
and  is  reported  as  a factor  of  safety  (F.S.).  The  factor 
of  safety  for  foundation  loadings  is  defined  as  the  ratio 
of  shear  strength  of  the  soil  to  the  shear  stress  required 
for  equilibrium  (Lambe  and  Whitman,  1969). 

F S = Shear  strength  of  soil  _ Equation  2-1 

Shear  stress  e>  equilibrium 

Therefore  when  the  shear  stress  equals  the  shear  strength 
(F.S.  = 1.0),  a gross  shear  failure  occurs  resulting  in 
excessive  deformations. 

The  use  of  this  limiting  equilibrium  definition  for 
the  factor  of  safety  requires  prediction  of  the  most  critical 
failure  surface  and  the  shear  strength  along  this  surface. 
There  are  more  than  twenty-five  methods  for  determining  the 
factor  of  safety  for  stability  analysis  (Whitman  and  Bailey, 
1967).  These  methods  vary  according  to:  the  assumed  shape 
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of  the  failure  surface  (circular,  log  spiral,  sliding  wedge 
or  general  surface);  the  method  of  determining  the  normal 
stress  along  this  failure  surface  and  thereby  determining 
the  shear  strength;  and  by  the  type  of  calculations  involved 
(arthimetic,  graphical  or  computer  solutions).  The  most 
popular  method  is  the  procedure  of  slices  as  first  pro'posed 
by  Fellenius  (1936),  which  is  also  known  as  the  Swedish  or 
Ordinary  Method  of  Slices.  This  method  has  been  adapted 
and  modified  to  varying  degrees  by  many  researchers,  e.g. 
Bishop  (1955),  Simplified  Bishop  (1955) » Janbu’s  Generalized 
Procedure  of  Slices  (1957),  Morgens tern-Price  (1965)  and 
Spencer's  method  of  slices  (196?) . These  methods  of  slices 
all  assume  some  shape  for  the  failure  surface,  then  divide 
the  soil  into  vertical  slices  and  apply  moment  and/or  force 
equilibrium  to  calculate  the  normal  forces  along  the  failure 
surface.  Various  assumptions  must  be  made  to  reduce  the 
unknowns  and  make  the  equations  statically  determinate. 

These  assumptions  involve  side  force  magnitude,  direction 
and  location  of  its  line  of  thrust  on  the  side  of  each  slice. 
Wright  (1969)  gives  an  excellent  summary  of  twenty-one 
different  limiting  equilibrium  methods  for  the  factor  of 
safety  and  compares  their  results  (See  Table  2-1). 

Generally  the  methods  agree  within  about  ten  percent  except 
for  the  Swedish  method  in  granular  soils  and  the  sliding 
block  method  in  cohesive  soils,  both  of  which  may  give 
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factors  of  safety  significantly  too  low. 

Sources  of  error  in  calculating  a factor  of  safety 
include  (Barboteu,  1972  and  Gilbert,  1974) » the  random 
testing  of  a variable  soil,  modeling  the  soil  as  layers  of 
homogeneous  material,  the  selection  of  the  most  critical 
failure  surface,  the  effects  of  end  restraint  in  actual 
three  dimensional  situations  (Baligh  and  Azzouz,  1975)*  and 
testing  procedures  used  to  evaluate  the  strength  parameters. 
The  testing  procedures  may  introduce  error  due  to  sample 
disturbance,  triaxial  versus  plane  strain  testing,  strain 
rates, soil  anisotropy  and  determination  of  the  changes  of 
the  in  situ  principle  stresses  during  loading.  Therefore, 
the  input  parameters  of  soil  strength  and  geometry  become 
more  important  sources  of  error  than  the  choice  of  method. 

For  example,  Johnson  (1974)  concludes » 

"The  use  of  total  stress  versus  effective  stress 
analyses  and  the  various  ways  in  which  design  strengths 
oan  oe  selected  produce  a wide  range  of  safety  factors 
and  are  more  important  than  the  method  for  analyzing 
stability. " 

The  method  used  for  stability  analysis  of  an  embankment 
is  up  to  the  user,  based  on  one's  own  experience;  choice  of 
arithmetic,  graphical  or  computer  solutions  and  the 
availability  of  computer  programs.  Even  the  most 
sophisticated  computer  program  cannot  insure  accuracy 
greater  than  the  input  parameters  nor  should  the  computer 


results  be  used  without  being  checked  for  reasonable 
solutions  (Little  and  Price,  1958). 

In  this  thesis  the  methods  used  are  the  Simplified 
Bishop  and  the  STAB3D  methods  of  slices.  The  Simplified  or 
Modified  Bishop  method  of  slices  as  programed  on  LEASE  II 
(Limiting  Equilibrium  Analysis  in  Soils  Engineering)  ("Dawson, 
1972),  a subsystem  of  ICES  (Integrated  Civil  Engineering 
System) , is  a reasonably  accurate  method  of  stability  analysis 
and  well  suited  to  the  circular  arc  failure  analyzed  in  the 
Portsmouth,  New  Hampshire  Experimental  Test  Section  in 
Appendix  A.  Although  this  method  is  not  the  most  rigorous, 
its  error  is  less  than  eight  percent  when  compared  to  any 
more  rigorous  method  (Wright,  Kulhawy  and  Duncan,  1969) . 

STAB3D  is  a method  recently  developed  at  M.I.T.  which  includes 
the  ability  to  account  for  three  dimensional  effects  and 
strength  anisotropy  (Azzouz,  1977).  This  method  was  used 
in  Appendix  B for  the  analysis  of  a hypothetical  embankment 
constructed  on  a clay  foundation  having  the  anisotropic 
strength  properties  of  Boston  Blue  Clay. 


EQUILIBRIUM  CONDITIONS  SATISFIED  BY  VARIOUS  METHODS  OF  STABILITY  ANALYSIS 

(Wright  1969) 
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TABLE  2-1 


This  condition  of  equilibrium  is  implicitly  satisfied  as  a result  of  the 
direct  consideration  of  other  equilibrium  conditions. 

The  original  presentation  of  this  method  was  for  a circular  arc  only. 


CHAPTER  3 


TYPES  OF  ANALYSES 

The  'type'  of  stability  analysis  distinguishes  between 
the  choice  of  expressions  for  calculating  the  shear  strength 
of  the  soil.  The  total  stress  analysis  (TSA)  inputs  the 
shear  strength  directly  as  it  is  assumed  to  be  independent 
of  the  total  normal  stress  acting  on  the  failure  surface. 
However,  the  strength  may  vary  with  the  inclination  of  the 
failure  surface,  the  stress  history  of  the  clay  and  the 
level  of  strain.  The  effective  stress  analysis  (ESA)  inputs 
the  drained  strength  parameters  of  the  clay  (c  and  , the 
pore  pressure  (u)  and  the  total  stress  (tf")  to  calculate  the 
shear  strength  using  the  expression  in  Equation  3-1 • 

7^f  = c + (<r  - u)  tan  Equation  3*1 

Most  limiting  equilibrium  methods  of  stability  analysis 
and  all  methods  of  slices  may  be  used  with  either  type  of 
analysis.  In  general,  the  two  types  of  analyses  do  not 
agree  and  result  in  different  numerical  values  for  the  factor 
of  safety.  The  following  sections  will  discuss  the  definition 
of  strength,  stress  paths,  and  means  for  measuring  the  soil 
parameters  for  each  types  of  analysis.  Chapter  5 will 
compare  the  factors  of  safety  computed  by  both  types  of 
analyses  for  several  classes  of  stability  problems. 
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3.1  TOTAL  STRESS  ANALYSIS 


The  use  of  the  TSA  in  this  thesis  is  limited  to  cases 
where  the  clay  is  saturated.  For  soft  foundation  clays  upon 
which  embankments  are  constructed,  this  is  generally  true 
and  not  a serious  limitation  for  the  TSA. 

The  TSA  used  in  this  thesis  is  not  the  same  as  a 0=0 
analysis.  The  theoretical  basis  for  the  0=0  analysis  is 
given  by  Skempton  ( 1 948 ) and  shown  in  Figure  3-1  for  some 
typical  triaxial  test  results.  Recognizing  that  water  is 
nearly  incompressable  compared  to  the  soil  skeleton,  any 
change  in  total  stress  on  a saturated,  undrained  sample 
during  shear  results  in  an  equal  change  in  the  pore  pressure. 
Since  an  change  in  the  total  stress  is  carried  by  a like 
change  in  the  pore  fluid,  no  additional  stress  is  transferred 
to  the  soil  and  the  effective  stress  remains  constant.  The 
0=0  analysis  uses  the  maximum  shear  from  an  undrained  test  to 
determine  the  shear  strength  of  the  soil.  The  failure  plane 
corresponding  to  the  maximum  shear  is  always  orientated  at 
an  angle  of  45°.  Bishop  and  Bjerrum  (I960)  demonstrated 
for  a simple  example  excavation  that  the  0=0  analysis  gave 
the  same  results  as  an  effective  stress  analysis  for  the 
stability  but  different  locations  for  the  failure  surface. 

The  0=0  analysis  will  always  give  a failure  surface  at  4-5° 
while  the  actual  failure  occurs  at  an  angle  of  45°  + 0/2. 
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Because  the  0=0  analysis  gives  the  maximum  shear  and  a 
failure  surface  always  inclined  at  45°  to  the  principal  stress 
it  is  not  used  in  this  thesis.  Instead  the  shear  stress  on 
the  failure  surface  at  failure  iTff)  and  the  actual  failure 
surface  inclined  at  45°  + 0/2  to  the  principal  stress  are 
used  for  all  total  stress  analyses.  The  strength  of  the 
soil  is  still  independent  of  the  changes  in  total  stress  but 
the  strength  and  inclination  of  the  failure  surface  are  more 
accurately  represented. 

3.1-1  UNDRAINED  SOIL  BEHAVIOR 

When  determining  the  undrained  strength  of  a soil 
several  characteristics  of  soil  behavior  are  important  to 
consider.  These  include:  stress  history,  normalized  soil 
behavior,  anisotropy  and  strain  compatability.  Ladd  et  al. 
(1977)  present  an  excellent  summary  of  these  soil  behavior 
characteristics.  The  stress  history  of  the  soil  is  important 
because  of  the  effects  which  the  overconsolidation  ratio 
(OCR)  and  lateral  stress  ratio  (KQ)  have  upon  the  undrained 
strength.  Many  clay  soils  exhibit  similar  behavior  when 
compared  to  other  tests  on  the  same  soil  at  the  same  OCR. 

By  dividing  the  test  results  by  the  effective  consolidation 
stress  (S’vc)  they  may  often  be  represented  by  a single 
normalized  curve  which  makes  presenting  and  evaluating  the 
clay  behavior  easier.  The  use  of  these  normalized  curves 
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of  soil  behavior  led  to  the  developement  of  the  SHANSEP 
(Stress  History  And  Normalized  Soil  Engineering  Properties) 
(Ladd  and  Foott,  19?4)  for  determining  the  undrained  strength. 
The  SHANSEP  method  uses  the  stress  history  of  the  soil  i.e. 
the  existing  effective  consolidation  stress  and  the  maximum 
past  pressure,  to  determine  the  OCR  and  from  the  plot" of 
normalized  undrained  strength  vs  OCR  to  select  the  appropriate 
normalized  strength.  This  normalized  strength  is  then 
multiplied  by  the  effective  consolidation  stress  to  obtain 
the  undrained  strength.  Appendix  B demonstrates  the 
application  of  the  SHANSEP  method  for  a hypothetical  embankment 
constructed  on  a soil  having  the  properties  of  Boston  Blue 
Clay. 

Soil  anisotropy  may  be  considered  to  consist  of  two 

components:  inherent  and  stress  induced  anisotropy.  The 

inherent  anisotropy  is  due  to  soil  fabric  or  orientation  of 

particles  which  result  from  deposition  of  the  soil.  Varved 

clays  are  an  example  of  a soil  which  even  macroscopically 

have  variation  in  the  soil  properties  with  orientation,  but 

almost  all  soils  have  some  degree  of  inherent  anisotropy. 

Stress  induced  anisotropy  is  present  whenever  the  lateral 

stress  ratio  (K  ) is  other  than  one.  This  stress  induced 
o 

anisotropy  is  not  measured  correctly  unless  the  in  situ 
lateral  stress  ratio  is  duplicated  in  the  test.  Of  course 
it  is  the  combination  of  both  components  of  anisotropy 
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which  is  important  in  determining  the  undrained  strength 

and  is  often  represented  by  K , the  anisotropic  strength 

s 

ratio  (Kg  = (H)/ 'T'ff  ( V) ) . This  ratio  of  horizontal  to 

vertical  shear  strength  may  be  as  low  as  0.3  and  must  be 
considered  in  any  accurate  determination  of  undrained  strength. 

Just  as  the  strength  of  a clay  varies  with  changes  in 
failure  mode,  so  do  the  strains.  Because  the  peak  shear 
strengths  occur  at  different  levels  of  strain  for  different 
failure  modes  (i.e.  plane  strain  compression,  direct  simple 
shear  and  plane  strain  extension)  a stability  analysis  which 
considers  the  peak  strengths  to  occur  simultaneously  along 
the  failure  surface  will  overestimate  the  strength.  To 
correct  the  strengths,  their  values  might  be  averaged  at 
constant  strain  and  the  peak  strength  from  this  averaging 
used  in  the  analysis  rather  than  assuming  that  the  peak 
strengths  for  different  failure  modes  occur  simultaneously. 
Strain  compatibility  is  included  in  the  evaluation  of  the 
stability  for  the  embankment  in  Appendix  B. 

3.1-2  DETERMINATION  OF  UNDRAINED  STRENGTH 

There  are  three  techniques  for  determining  the  ur drained 
strength  of  a soil:  lab  tests,  in  situ  tests  and  empirical 
correlations.  For  lab  tests  there  are  two  general  approaches 
used  in  practice  to  determine  the  undrained  strength  for  a 
total  stress  analysis.  One  assumes  a unique  relationship 


between  the  water  content  at  failure  and  the  undrained 
strength  so  that  all  tests  are  conducted  at  the  natural 
water  content.  The  second  approach  assumes  a unique 
relationship  between  the  effective  consolidation  stress  and 
the  undrained  strength.  Tests  which  are  run  at  the  natural 
water  content  include i unconfined  compression,  "UU”  tr'iax.ial 
compression  (Bishop  and  Bjerrum,  I960)  Labvane  and  Torvane. 

I 

Tests  in  which  the  sample  is  reconsolidated  to  match  the 
effective  consolidation  stress  include i direct  shear,  direct 
simple  shear  and  "CU”  triaxial  compression  tests  (Ladd,  1971). 
By  conducting  laboratory  tests,  the  boundary  conditions  can 
be  controlled  ana  the  results  more  easily  interpreted  than 
in  situ  tests.  The  most  serious  drawback  to  lab  testing  is 
the  effect  of  sample  disturbance  which  can  reduce  the  measured 
undrained  strength  in  "UU"  tests  by  as  much  as  fifty  percent 
(Ladd  et  al.  1977).  The  "UU"  test  may  give  correct  solutions 
only  due  to  a fortuitous  cancelation  of  errors.  The 
reduction  in  strength  from  sample  disturbance  and  an  increase 
in  strength  from  the  rapid  shear  rate  and  use  of  the  strength 
due  to  vertical  loading  without  considering  anisotropy  may 
cancel  out  and  result  in  the  correct  strength.  By 
consolidating  the  soil  to  match  the  in  situ  lateral  stress 
ratio  (K0)  and  loading  beyond  the  maximum  past  pressure 
with  rebound  to  the  in  situ  CCR , Ladd  et  al.  (1977)  believe 
the  effects  of  sample  disturbance  are  minimized.  The  test 
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results  may  then  be  used  by  applying  normalized  strength 
properties  such  as  those  used  in  the  SHANSEP  method. 


In  situ  tests  offer  an  improvement  over  laboratory 
tests  because  of  their  generally  lower  costs  and  reduction 
in  the  disturbance  before  testing,  but  it  is  often  more 
difficult  to  determine  the  boundary  conditions  for  the  in 
situ  tests  and  the  results  may  be  more  difficult  to  interpret. 
The  field  vane  test  is  the  most  popular  field  test  for 
measuring  the  undrained  strength  but  its  failure  mode  is 
unique  and  its  strain  rate  much  higher  than  would  be 
experienced  under  an  embankment  loading.  The  Dutch  Cone 
test  is  becoming  increasingly  popular  but  its  results  are 
often  difficult  to  interpret  (Schmertmann,  1975). 


Empirical  correlations  are  often  used  in  soil  mechanics 


because  practice  usually  preceeds  the  theory  required  for  a 
thorough  analysis.  Most  empirical  correlations  are  based  on 
index  tests  such  as  Bjerrum's  (1972)  correction  factor  for 
field  vane  tests.  By  analyzing  the  circular  arc  failure  of 
sixteen  embankment  failures  based  on  the  field  vane  test,  he 
was  able  to  correlate  the  factor  of  safety  and  the  Plasticity 
Index  of  the  soil.  Although  Bjerrum  (1973)  explained  the 
correction  factor  by  a combination  of  strain  rate  effects 
and  anisotropy  for  different  plasticity  it  remains  as  an 


empirical  correlation.  The  Dutch  Cone  is  also  interpreted 


factors,  empirical  correlations  are  useful  as  checks  on  any 
lab  or  field  testing. 

3.2  EFFECTIVE  STRESS  ANALYSIS 

The  effective  stress  'type'  of  analysis  uses  the  drained 
strength  parameters  (c  and  ]5),  pore  pressure  (u)  and  the  total 
normal  stress  (<r)  to  express  the  shear  strength  of  the  soil 
(Equation  3.1).  This  expression  is  the  classical  Mohr 
Coulomb  failure  criteria  and  is  the  generally  accepted  way 
to  analyze  most  soil  engineering  stability  problems.  The 
ESA  is  often  preferred  over  the  TSA  because  it  is  known  that 
soil  behavior  is  controlled  by  the  effective  stress  and 
therefore  this  type  of  analysis  is  more  fundamental.  Also 
the  ESA  is  not  limited  to  saturated  soils  or  undrained 
loading  as  is  assumed  for  the  TSA.  Many  kinds  of  problems 
can  be  correctly  solved  by  the  ESA,  such  as  problems  in 
granular  soils,  unloading  in  cohesive  soils  and  the  stability 
of  natural  slopes.  The  "CD"  class  of  embankment  problems 
(defined  in  Section  4.2)  on  soft  clay  can  also  be  easily 
solved  using  the  ESA. 

To  use  the  ESA  the  failure  envelope  must  be  defined  by 
the  effective  strength  parameters  (c  and  0’)  and  the  changes 
in  the  total  stress  and  pore  pressure  must  be  known.  The 
changes  in  total  stress  are  usually  determined  by  applying 
one  of  the  methods  of  slices  to  compute  the  stress  normal  to 
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the  failure  surface.  The  evaluation  of  the  failure  envelope 
is  usually  made  by  conducting  shear  tests  on  soil  samples 
reconsolidated  in  the  laboratory.  The  tests  can  be  either 
CD,  consolidated-drained,  or  CU,  consolidated-undrained  with 
pore  pressure  measurements.  The  CD  tests  include  direct  box 
shear  and  triaxial  compression  tests.  The  CU  test  ard 
consolidated-undrained  test  in  which  the  pore  pressures  have 
been  measured  throughout  the  test  so  the  effective  stresses 
can  be  calculated.  Bishop  and  Bjerrum  (i960)  give  an 
excellent  review  of  all  triaxial  tests,  their  applications 
and  sources  of  error.  Ladd  (1971)  and  Ladd  et  al.  (1977) 
discuss  the  difficulties  and  sources  of  error  for  all  tests. 
Possible  sources  of  error  in  triaxial  tests  includes  the 
effects  of  sample  disturbance;  end  restraints;  pore  pressure 
equilization;  correction  factors  for  filter  strips,  membrane 
and  piston  friction;  back  pressure  for  saturation;  efficency 
of  filter  strips;  anisotropy;  and  test  interpretation.  The 
friction  angle  (^)measured  at  the  maximum  stress  (qf)  is 
lower  than  the  friction  angle  measured  at  maximum  obliquity 
(d' • Anisotropy  may  also  be  significant  especially  in 
lean  sensitive  clays  where  the  friction  angle  at  qf  from 
vertical  loading  may  be  as  much  as  ten  degrees  lower  than 
the  friction  angle  for  horizontal  loading  (Ladd  et  al.  1977). 
The  cohesion  intercept  (c)is  often  smaller  than  the  correction 
factor  for  filter  strips  although  for  loadings  on  soft  clay 
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the  cohesion  can  safely  be  taken  to  be  zero  for  most  soils. 

The  most  difficult  parameter  to  determine  for  an  ESA 
is  the  pore  pressure.  During  construction  piezometers  can 
be  installed  throughout  the  foundation  clay  and  measurements 
from  the  piezometers  used  to  create  a pore  pressure  field 
for  evaluation  of  the  effective  stresses.  This  technique 
can  be  used  to  relate  directly  the  piezometer  readings  to 
the  factor  of  safety  for  the  embankment.  This  approach  is 
attractive  for  many  engineers  but  it  does  not  give  any 
indication  of  changes  in  pore  pressure  due  to  shear.  This 
technique  of  using  piezometers  to  determine  the  pore  pressure 
beneath  the  embankment  does  not  help  the  engineer  designing 
the  embankment  and  only  serves  as  a tool  to  monitor  construc- 
tion. To  predict  pore  pressure  beneath  the  embankment  prior 
to  any  construction  is  extremely  difficult  and  even  with 
pore  pressure  models  such  as  YLIGHT  (Section  4.1)  it  is  not 
yet  possible  to  accurately  predict  the  pore  pressures 
generated  throughout  an  embankment  foundation. 

Figure  3-2  shows  a stress  path  to  failure  for  a "typical 
element"  beneath  an  embankment  on  soft  clay.  When  computing 
the  shear  strength  of  the  soil  using  the  expression  in 
Equation  3*1  for  the  ESA  one  obtains  the  shear  on  the 
failure  surface  at  a constant  effective  stress.  For 
normally  consolidated  clays,  this  strength,  TffESA,  is 
higher  than  the  strength  obtained  from  a TSA  because  it 
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does  not  include  the  increased  pore  pressures  generated 
during  shear  to  reach  the  failure  envelope.  As  the  state 
of  stress  on  the  "typical  element"  approaches  the  failure 
envelope  the  differences  between  the  two  types  of  analysis 
will  decrease.  At  failure  conditions  both  the  ESA  and  TSA 
should  give  the  same  value  for  the  strength. 

In  general,  the  strengths  computed  by  the  two  types  of 
analyses  are  not  the  same  and  for  the  case  of  normally 
consolidated  clay  the  ESA  gives  higher  strengths  than  does 
the  TSA.  Therefore,  the  factor  of  safety  for  stability 
computed  by  the  two  types  of  analyses  are  different  because 
the  values  of  strength  used  in  the  limiting  equilibrium 
definition  for  the  factor  of  safety  are  different.  The  ESA 
should  give  a factor  of  safety  higher  than  the  TSA  for  all 
embankments  on  soft  clay  except  at  failure  when  both  should 
give  a factor  of  safety  of  1.0.  The  two  types  of  analyses 
are  compared  in  Chapter  5 for  several  cases  of  embankment 
stability. 
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CLASSES  OF  STABILITY 

Any  construction  (loadings  or  excavations)  on  clay 
soils  results  in  changes  in  pore  pressure.  This  change  is 
often  assumed  to  occur  instantaneously  with  construction 
and  then  vary  with  time  until  the  pore  pressures  have 
reached  equilibrium  or  steady  state  conditions.  To  evaluate 
the  stability  of  construction  on  clay  it  is  convenient  to 
define  three  'classes'  for  stability  based  on  the  changes 
in  pore  pressure  due  to  drainage. 

The  first  class  is  the  case  when  no  drainage  occurs  and 
hence  there  has  been  no  dissipation  of  excess  pore  pressures 
during  loading  or  unloading.  It  is  frequently  assumed  that 
construction  proceeds  rapidly  enough  to  satisfy  this 
condition  and  hence  this  class  is  known  as  the  "UU"  = 
unconsolidated-undrained,  end-of-construction,  immediate  or 
short  term  stability. 

The  second  class  of  stability  is  when  the  clay 
foundation  has  reached  equilibrium  or  steady  state  conditions 
wherein  consolidation  is  complete  and  there  are  no  excess 
pore  pressures.  Failure  is  then  assumed  to  occur  slowly 
without  any  excess  pore  pressure  developed  during  shear. 

This  class  is  known  as  "CD"  = consolidated-drained  or  long 
term  stability. 
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The  third  class  is  intermediate  between  'end-of- 
construction'  and  'long  term’,  when  there  has  been  some 
degree  of  drainage  so  excess  pore  pressures  have  partially 
dissipated  but  not  yet  reached  equilibrium  conditions  and 
failure  occurs  so  rapidly  that  there  is  essentially  no 
drainage.  Alternately,  the  foundation  clay  may  be  fully 
consolidated  under  the  applied  loads,  but  failure  again 
occurs  so  rapidly  that  excess  pore  pressures  generated  by 
the  failure  do  not  have  time  to  dissipate.  This  class  is 
known  as  "CU"  = consolidated-undrained,  intermediate  or 
partial  drainage  stability. 

In  theory,  each  of  the  three  'classes'  of  stability 
could  be  analyzed  by  either  'type'  of  stability  analyses, 
i.e.  TSA  or  ESA.  However  in  practice  one  type  of  analysis 
may  be  much  better  suited  to  handle  a particular  class  of 
problem.  The  following  sections  will  review  each  class  of 
stability  by  both  the  total  stress  and  the  effective  stress 
analyses  and  discuss  the  advantages  and  disadvantages  of 
each  for  embankments  on  soft  clay. 

4.1  "UU"  END-OF-CONSTRUCTI ON 

The  "UU"  class  of  stability  problems  is  the  most 
critical  for  embankments  on  soft  clays  where  the  applied 
stresses  are  such  that  the  foundation  clay  is  loaded  beyond 
the  in  situ  maximum  past  pressure  of  the  soil.  (Note  that 
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for  excavations  and  natural  slopes  the  "CD"  class  is  most 
critical).  An  embankment  on  a deep  truely  normally 
consolidated  clay  is  usually  constructed  much  faster  than 
excess  pore  pressures  can  dissipate  so  there  is  little  or 
no  error  in  assuming  undrained  conditions.  However, 
embankments  built  slowly,  over  a relatively  thin  layef  of 
clay  or  on  an  overconsolidated  clay  would  not  be  an  undrained 
loading  condition  and  must  be  treated  as  a "CU"  class 
stability  problem.  In  any  case,  if  the  clay  is  saturated 
and  then  loaded  quickly,  the  assumptions  for  a conventional 
total  stress  analysis  have  been  met  and  the  "UU"  class  of 
stability  has  usually  been  solved  using  the  total  stress 
analysis.  The  TSA  for  "UU"  class  embankment  problems  is 
recommended  by  Bishop  and  Bjerrum  (i960),  Ladd  (1971)  and 
Janbu  (1977). 

"However,  the  main  reason  for  using  a total  stress 
approach  for  short-term  undrained  conditions  in 
clay  is  to  avoid  predicting  the  pore  pressures  and 
hence  the  effective  stresses.  This  may  at  times  be 
justified  for  embankments  and  footings  when  the  load 
on  the  clay  exceeds  the  preconsolidation  load." 

(Janbu,  1977) 

The  only  reason  for  any  further  discussion  concerning 
the  TSA  for  the  "UU"  class  of  stability  is  that  the  TSA 
doesn't  always  give  the  correct  factor  of  safety  and 
embankment  failures  have  resulted.  The  problem  is  not 
with  the  method  of  analysis  or  the  total  stress  approach 


to  the  problem  but  in  the  value  of  the  undrained  strength 
used  in  the  calculations.  For  example,  Bjerrum  (1972) 
concludes: 

"Obviously,  there  is  nothing  wrong  with  the  general 
principle  of  the  presently  used  methods  of  computing 
the  stability.  When  they  do  not  work,  it  is  the  shear 
strength  introduced  in  the  computation  which  is  * 
incorrect. " 

To  correct  input  strength  based  on  field  vane  measurements 
Bjerrum  (1972)  recommended  an  empirical  correction  for  field 
vane  tests.  As  stated  above,  it  does  not  imply  any  correction 
to  the  total  stress  analysis  but  only  to  the  measurement  of 
the  undrained  shear  strength.  The  difficulties  of  determining 
the  undrained  shear  strength  from  other  tests  were  discussed 
in  Section  3*1 • 

The  effective  stress  analysis  is  not  often  used  for 
the  "UU"  class  of  stability  problems  because  of  the  difficulty 
in  predicting  pore  pressures.  Also  the  "UU"  class  implies 
undrained  conditions  to  failure,  whereas  the  effective 
stress  analysis  uses  a drained  stress  oath  to  determine  the 
strength  at  failure.  The  ESA  does  not  include  any  change 
in  the  effective  stress  due  to  shear  to  failure  as  the 
existing  effective  stress  used  in  the  equation  for  failure 
strength  (Equation  3*1)  results  in  a failure  condition  at 
the  same  effective  stress.  Schmertmann  (1975)  recommends 
the  use  of  the  ESA  exclusively  because  of  the  insight  into 


the  basic  soil  behavior  and  because  the  empirical  correction 
factors  used  in  obtaining  the  undrained  strength  for  the  TSA 
e.g.  Bjerrum's  field  vane  correction,  have  no  theoretical 
justification.  Perhaps  when  testing  technology  advances  so 
that  the  undrained  shear  strength  is  more  accurately  measured 
his  argument  against  using  a TSA  will  be  eliminated. 

For  embankments  constructed  on  lightly  overconsolidated 
clays  (0CR<2.5)t  Lero  ieil  et  al.  (1978)  have  shown  that 
the  foundation  clay  acts  as  a partially  drained  material, 
due  to  the  high  coefficient  of  consolidation  cv,  until  the 
preconsolidation  stress  is  reached  at  which  time  the  clay 
becomes  normally  consolidated  and  acts  as  an  undrained 
material.  The  effects  of  drainage  during  the  recompression 
of  the  overconsolidated  clay  would  usually  result  in  a 
slight  increase  in  strength  and  the  analysis  based  on  the 
initial  in  situ  undrained  strength  would  be  conservative. 
However,  with  highly  structured,  sensitive  clays  (such  as 
the  Lake  Champlain  clays  in  Canada)  recompression  may 
actually  result  in  a decrease  in  undrained  strength  according 
to  Tavenas  et  al.  (1978). 

For  lighly  overconsolidated  clays  Leroueil  et  al.  (1978) 
recommend  prediction  of  pore  pressures  based  upon  the  YLIGHT 
(Yield  Locus  Influenced  by  Geological  History  and  Time) 
model  of  clay  behavior  (Figure  4-1).  The  YLIGHT  model 
predicts  a rate  of  pore  pressure  generation  of  = 0.6, 
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where  B is  the  ratio  of  the  change  in  pore  pressure  to  the 
change  in  total  vertical  stress,  until  the  preconsolidation 
load  (P  ) is  reached  at  the  critical  height  (Hcr) . The 
foundation  clay  becomes  normally  consolidated  under  the 
embankment  load  at  the  critical  height  and  further  loading 
generates  pore  pressure  at  the  rate  of  B2  = 1.0  in  a true 
undrained  loading  condition.  At  conditions  of  local  yielding 
with  a strain  softening  material  the  pore  pressures  are 
generated  at  a rate  in  excess  of  one.  Though  this  model  is 
very  helpful  in  depicting  pore  pressure  generation  in  lightly 
overconsolidated  clays  it  is  only  valid  for  the  center  line 
of  the  embankment  and  cannot  be  used  to  predict  pore  pressure 
throughout  the  entire  foundation  which  is  necessary  for  an 
effective  stress  analysis. 

4.2  "CD"  LONG  TERM 

The  long  term  class  of  stability  assumes  complete 
drainage  with  no  excess  pore  pressures  due  to  construction 
or  shear  to  failure.  For  embankments  on  soft  clay  the  long 
term  is  never  the  critical  class  for  stability  analysis. 

With  consolidation  under  the  embankment  load  the  effective 
stress  in  the  clay  increases,  resulting  in  an  increased 
shear  strength  so  that  the  factor  of  safety  is  also  increased. 
Often  the  factor  of  safety  for  this  class  of  stability  is 
calculated  only  to  insure  that  it  is  above  some  preset  value 


arbitrarily  selected  prior  to  design.  Since  this  class  of 
stability  assumes  complete  consolidation  and  no  change  in 
effective  stresses  during  shear  to  failure  i.e.  they  remain 
equal  to  the  consolidation  stresses,  there  are  no  excess 
pore  pressures.  The  equilibrium  pore  pressures  can  therfore 
be  determined  from  the  water  table  or  if  steady  state  'flow 
exists  from  a flow  net.  Figure  4-2  shows  the  stress  path 
for  a point  on  the  failure  surface  of  a "typical  element" 
beneath  an  embankment!  0-1  represents  the  loading  due  to 
construction  and  1-2  the  consolidation. 

The  effective  stress  analysis  is  used  almost  exclusively 
for  this  class  of  stability  problems  because  it  is  easily 
applied.  Since  there  are  no  excess  pore  pressures  the  ESA 
is  the  same  as  assuming  a drained  loading  from  2-ESA2,  Figure 
4-2.  The  use  of  the  ESA  for  the  "CD"  class  of  stability  is 
easy  to  apply  but  of  little  practical  significance  for 
embankments  (except  for  those  constructed  in  stages,  see 
Section  4.3)  because  the  higher  strengths  due  to  consolidation 
make  the  ’ end-of-construction'  more  critical  for  design. 

If  a fully  consolidated  embankment  was  quickly  loaded 
i.e.  additional  construction  due  to  enlargement  of  an 
existing  embankment  or  an  extraordinarily  heavy  transient 
loading,  the  additional  load  would  produce  an  undrained 
loading,  2-TSA2  in  Figure  4-2.  A total  stress  analysis 
would  now  be  appropriate  for  thsse  conditions  of  additional 


load.  The  question  is  whether  a rapid  failure  or  the  very 
slow  drained  failure  assumed  in  the  ESA  better  represents 
the  actual  long  term  stability  of  an  embankment.  Treating 
the  long  term  stability  with  undrained  loading  to  failure 
is  the  same  as  a "CU"  class  stability  problem  with  complete 
consolidation.  This  comparison  will  be  discussed  further 
in  Section  5»3» 

4.3  "CU"  INTERMEDIATE 

The  ''CU"  class  of  stability  can  apply  to  two  conditions 
of  embankment  construction.  One  of  these  is  when  some 
drainage  has  occurred  but  the  consolidation  is  not  yet 
complete.  If  the  drainage  under  the  embankment  were  complete 
the  case  becomes  the  same  as  the  long  term  class  assuming 
that  no  excess  pore  pressures  are  developed  during  actual 
failure,  and  if  the  clay  is  completely  undrained  it  is  the 
same  as  a "UU"  class  stability  problem.  This  would  be  the 
same  as  determining  the  variation  in  the  factor  of  safety 
with  the  drainage  time,  with  the  "UU"  class  for  the  initial 
factor  of  safety  when  there  is  no  drainage,  the  long  term 
class  for  completely  consolidated  and  the  "CU"  class  for 
any  time  in  between  these  extremes. 

The  second  condition  for  the  "CU"  class  is  the  staged 
construction  of  embankments  wherein  construction  is 
interrupted  to  allow  additional  drainage  before  continuing 
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construction.  This  interruption  in  work  could  be  due  to 
seasonal  weather  changes,  work  strikes,  modification  to  an 
existing  embankment  or  an  intentional  delay  to  allow  the 
clay  beneath  the  embankment  time  to  consolidate  and  gain 
strength  so  the  embankment  may  be  built  higher  than  could 
be  allowed  in  a single  stage. 

A total  stress  analysis  for  the  "CU"  class  of  stability 
implies  an  undrained  failure  starting  from  the  existing 
effective  consolidation  stress.  As  the  c^ay  consolidates 
under  the  embankment  load  the  excess  pore  pressures  dissipate 
resulting  in  increased  effective  stresses  and  increased 
undrained  strength.  The  two  major  difficulties  in  this  type 
of  analysis  are  estimating  or  measuring  the  existing 
effectives  stresses  beneath  the  embankment  and  relating  the 
undrained  strength  to  these  consolidation  stresses.  One 
way  to  obtain  the  existing  effective  consoli dation  stresses 
in  the  clay  foundation  is  to  evaluate  the  preconsolidation 
stress  history  of  the  clay,  the  changes  in  total  stress  and 
pore  pressure  due  to  loading  and  the  change  in  excess  pore 
pressure  due  to  consolidation  drainage.  Once  the  existing 
stresses  are  determined  the  undrained  strength  could  be 
determined  by  applying  normalized  soil  properties. 
Alternately,  one  might  directly  measure  the  existing 
undrained  strength  by  in  situ  testing  or  undisturbed 
sampling  and  lab  shear  tests. 
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Turnbull  and  Hvorslev  (1966)  recommend  the  use  of  the 
total  stress  analysis  for  the  "CU"  class  of  stability 
problems,  and  state  1 

’’The  appropriate  shear  strength  is,  therefore,  the 
undrained  strength  of  the  soil  after  it  has  been 
consolidated  to  a condition  corresponding  to  the 
estimated  effective  stresses  in  the  field.” 

Ladd  (1971)  recommends  the  use  of  the  TSA  with  SHANSEP 
method  for  this  class  of  problems.  Tavenas  et  al.  (1978) 
used  field  vane  tests  after  partial  consolidation  to  measure 
the  increased  strength  for  the  TSA  in  this  class  of  stability. 

The  effective  stress  analysis  for  the  "CU"  class  of 
stability  problems  requires  knowledge  of  the  existing 
effective  consolidation  stresses  just  as  in  the  TSA.  Failure 
is  then  assumed  to  occur  without  generation  of  additional 
excess  pore  pressures.  It  is  inconsistent  to  consider  the 
existing  excess  pore  pressures  to  remain  constant  while 
allowing  the  excess  pore  pressures  generated  during  shear 
to  dissipate  as  is  assumed  for  all  ESA. 

To  determine  the  existing  effective  consolidation 
stresses  beneath  the  embankment  one  could  subtract  the  loss 
of  pore  pressure  due  to  consolidation  drainage  from  the 
pore  pressure  generated  during  loading  as  mentioned  in  the 
TSA  above.  Another  way  is  to  measure  pore  pressures 
throughout  the  foundation  with  piezometers  and  use  the 
measured  pore  pressures  to  determine  the  existing  effective 
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consolidation  stresses.  The  effective  stress  analysis  for 
stability  of  an  embankment  can  then  be  directly  related  to 
the  pore  pressures  measured  from  piezometers  installed 
throughout  the  foundation  clay.  This  ability  to  monitor 
pore  pressures  and  calculate  the  effective  stress  factor  of 
safety  for  any  instant  in  time  makes  the  ESA  attractive  to 
many  engineers.  However,  the  factor  of  safety  calculated 
does  not  include  changes  in  pore  pressure  resulting  from 
shear  due  to  actual  failure.  Since  pore  pressures  may 
increase  rapidly  near  failure  loads,  a larger  margin  of 
error  must  be  allowed  when  using  the  ESA.  For  example 
Tavenas  et  al.  (1978)  state: 

"When  some  points  of  the  foundation  reach  local  failure 
. . . increased  pore  pressures  are  generated  and  the 
available  strength  decreases  ...  as  a result  the 
factor  of  safety  will  decrease  rapidly  toward  1.0  at 
complete  failure." 

Figure  4-3  depicts  the  generalized  behavior  predicted  by 
Tavenas  et  al.  (1978)  with  the  ESA  factor  of  safety  decreasing 
rapidly  toward  1.0  after  reaching  local  yielding  near  a factor 
of  safety  of  1.3.  This  behavior  is  discussed  further  in 
Section  5*2. 

The  ESA  for  the  "CU"  class  of  stability  problems  is 
recommended  by  Schmertmann  (1975) » Janbu  (1977)  and  Tavenas 
et  al.  (1978) . 
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B = Rate  of  pore  pressure  generation  = AU /t^c \ 

H = Height  of  embankment. 

Hcr  = Critical  height  when  foundation  clay  becomes 
normally  consolidated. 

I.  = Influence  factor  for  loading. 

PQ  = Initial  vertical  effective  stress. 

Pc  = Maximum  past  pressure. 

U = Pore  pressure. 

= Total  unit  weight  of  embankment  soil. 
a ^ = Total  vertical  stress. 


YLIGHT  MODEL  OF  PORE  PRESSURE  PREDICTION 

(Leroueil  et  al.  19?3) 

FIGURE  4-1 
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CHAPTER  5 


COMPARISON  OF  TOTAL  vs  EFFECTIVE  STRESS  ANALYSES 

The  true  test  of  either  'type'  of  stability  analysis 
i.e.  TSA  or  ESA,  for  embankments  on  soft  clay  is  its  ability 
to  correctly  predict  a failure  condition.  If  both  types  of 
analyses  give  an  equal  degree  of  accuracy  for  the  prediction 
of  failure  conditions  then  the  choice  of  which  type  of 
analysis  to  use  is  up  to  the  engineer.  However  if  one 
approac'  can  be  demonstrated  to  be  significantly  more 
accurate,  then  that  approach  should  be  used  over  the  other 
type  of  analysis.  In  an  effort  to  make  this  comparison 
between  TSA  and  ESA  for  embankments  on  soft  clay,  a review 
of  the  literature  was  necessary  to  compile  existing  data. 
Despite  numerous  embankment  stability  case  studies  reported 
in  the  literature,  there  are  very  few  authors  who  have 
calculated  stability  by  both  the  TSA  and  ESA  for  the  same 
embankment.  This  may  be  due  to  a bias  toward  one  type  of 
analysis  by  the  authors  and/or  due  to  difficulties  in  using 
one  type  of  analysis  for  a particular  class  of  stability. 

The  following  sections  will  compare  the  TSA  and  ESA  of 
embankments  on  soft  clay  for  three  conditions:  actual  failures, 
variation  in  the  factor  of  safety  with  embankment  height  and 
the  long  term  stability. 


5.1  CASES  OF  EMBANKMENT  FAILURE 

At  failure  conditions  the  value  of  the  limiting 
equilibrium  factor  of  safety,  as  defined  in  Chapter  2,  is 
1.0.  By  comparing  the  results  of  both  types  of  analyses 
for  failure  conditions,  the  values  of  the  factor  of  safety 
can  be  compared  to  the  known  value  of  1.0  and  not  just  to 
the  other.  For  embankment  conditions  other  than  failure, 
the  actual  factor  of  safety  is  unknown  and  a comparison  of 
the  factor  of  safety  from  a TSA  and  an  ESA  would  be  only 
relative . 

The  ' end-of-construc tion'  class  of  stability  is  the 
most  critical  class  for  embankments  on  soft  clay  and  these 
problems  have  usually  been  treated  using  a TSA  (Section  ^.1). 

A review  of  embankment  failures  reported  in  the  literature  has 
produced  eight  cases  which  were  analyzed  by  both  types  of 
analyses.  These  eight  embjnkment  failures  were  reported  by 
a total  of  twelve  researchers  and  are  listed  in  Table  5-1 • 

Two  of  the  embankments,  Portsmouth  and  1-95 » were  test 
sections  intentionally  constructed  to  failure.  The 
Atchafalaya  levee  didn't  fail  by  a rapid  rotation  or  sliding 
motion  but  the  undrained  deformations  were  large  and  tension 
cracks  appeared  in  the  fill  so  that  both  researchers 
concluded  that  the  embankment  was  at  failure.  The  other 
five  embankments  experienced  accidental  failures  during 
construction.  Considering  the  large  number  of  embankment 
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failures  reported  in  the  literature  this  small  number  of 
embankments  analyzed  with  both  types  of  analysis  is  very 
small  and  perhaps  indicates  how  seldom  the  ESA  is  used  for 
the  "UU"  class  stability  problem. 

The  TSA  for  these  embankment  failures  are  divided  into 
six  groups  based  upon  the  technique  used  to  determine  'the 
undrained  strength.  These  groups  are:  field  vane,  corrected 
field  vane  (Bjerrum,  1972),  SHANSEP  (Ladd  and  Foott,  197^)« 
unconfined  compression,  unconsolidated-undrained  triaxial 
compression  tests  and  an  average  of  the  unconfined  and  "UU" 
triaxial  compression  tests.  The  ESA  is  divided  into  two 
groups  based  upon  the  method  of  determining  the  pore  pressure, 
by  measuring  piezometer  levels  or  predictions  based  on  the 
loading  conditions.  All  the  results  are  given  in  Table  5-1 
and  shown  graphically  in  Figure  5-1* 

Both  types  of  analyses  show  a wide  fluctuation  in  the 
computed  factors  of  safety  and  considering  all  data  equally 
neither  type  seems  to  be  very  accurate.  However  some  of 
the  techniques  used  to  measure  the  undrained  strength  are 
known  to  introduce  a large  source  of  error  into  the  calcula- 
tions. If  these  techniques  and  their  resulting  factors  of 
safety  are  eliminated  from  the  comparison  the  accuracy  of 
the  TSA  is  vastly  improved  and  certainly  superior  to  the 
ESA.  The  field  vane  strengths  have  been  shown  by  Bjerrum 
(1972)  to  overestimate  the  shear  strength  of  plastic  soil 


and  the  field  vane  strengths  should  not  be  used  without 

applying  Bjerrum's  empirical  correction  factor.  The  corrected 

field  vane  values  of  the  factor  of  safety  are  quite  good 

except  for  the  Portland  Grganic  Clay.  This  clay  had  wood 

pieces  and  shells  throughout  the  layer  and  the  field  vane 

results,  even  with  the  correction  factor,  would  be  expected 

to  give  unacceptably  high  results.  The  unconfined  compression 

tests  can  also  introduce  large  amounts  of  error  due  to 

sample  disturbance  and  should  not  be  considered  accurate  in 

any  case.  With  all  of  these  TSA  test  results  eliminated 

from  consideration  and  only  the  lower  value  of  any  range 

given  by  a autnor  used  in  the  comparison,  the  range  of  the 

TSA  based  on  SHANSEP  is  1.08  to  0.82  with  an  average  of  0.93, 

the  TSA  based  on  corrected  field  vane  tests  is  1.30  to  0.99 

with  an  average  of  1.0?  and  the  ESA  ranges  from  1.95  to  0.96 

with  an  average  of  1.2?.  These  results  from  actual  embankment 

failures  plotted  in  Figure  5-2  support  the  use  of  TSA  for  the  j 

"UU"  class  of  stability  for  embankments  on  soft  clay.  I 


5.2  VARIATION  WITH  EMBANKMENT  HEIGHT 

A comparison  of  the  factor  of  safety  calculated  by  both 
types  of  analyses  at  any  condition  other  than  failure  is 
less  meaningful  because  the  actual  factor  of  safety  is 
unknown.  Nevertheless  the  variation  of  the  factor  of  safety 
with  changes  in  the  embankment  height  should  show  a trend 


< 


I 


of  decreasing  factor  of  safety  with  increasing  height.  At 
the  failure  condition  the  f .ctor  of  safety  for  both  types 
of  analyses  should  be  1.0. 


To  determine  the  variation  of  the  factor  of  safety 
with  embankment  height  of  an  embankment  constructed  to 
failure,  the  Portsmouth,  New  Hampshire  Experimental  Te"st 
Section  (Ladd,  1972)  was  selected.  The  Portsmouth  embankment 
was  constructed  on  a soft  sensitive  clay  and  heavily 
instrumented  so  both  a TSA  and  an  ESA  could  be  conducted. 

The  soil  properties,  embankment  geometry  and  stability 
analyses  are  given  in  Appendix  A and  will  not  be  repeated 
here.  The  TSA  was  based  on  the  average  SHANSEP  strength 
values  computed  for  CKQU  tests  for  shear  and  oedometer  tests 
for  stress  history.  The  ESA  was  based  on  measured  pore 
pressures  and  the  drained  strength  parameters,  = 21°  and 
c = 0,  which  were  the  lower  bound  from  CIU  tests  on  the  soil. 
The  values  of  the  factor  of  safety  with  changes  in  the 
embankment  height  are  listed  in  Tables  A-4  and  A-5  of 
Appendix  A and  shown  graphically  in  Figure  5~3*  The  factor 
of  safety  from  an  effective  stress  analysis  is  always  higher 
than  that  from  a total  stress  analysis  but  their  values 
approach  each  other  near  the  failure  height. 

In  addition  to  the  Portsmouth  embankment,  six  other 
embankments  are  listed  in  Table  5~2  where  the  variation  of 


both  the  TSA  and  ESA  factor  of  safety  with  embankment  height 


were  reported.  The  TSA  for  the  four  embankments  reported 
by  Tavenas  et  al.  (1978)  were  computed  using  uncorrected 
field  vane  tests  after  each  increment  of  loading  for  these 
stage  constructed  embankments.  The  TSA  for  the  New  Liskeard 
varved  clay  embankment  were  based  on  average  unconfined  and 
unconsolidated-undrained  triaxial  compression  tests.  'The 
other  embankments  listed  in  Table  5-2  used  the  SHANSEP 
technique  for  calculating  the  TSA  factor  of  safety.  All  of 
the  embankments  had  piezometers  to  measure  the  pore  pressure 
which  were  used  to  determine  the  ESA  factor  of  safety.  The 
results  for  these  seven  embankments  are  shown  in  Figure  5-4 
which  plots  the  TSA  factor  of  safety  versus  the  ESA  factor 
of  safety  for  different  elevations  of  the  embankments.  In 
every  case  the  ESA  gives  a higher  value  for  the  factor  of 
safety  than  the  TSA  for  the  same  height  of  embankment. 

With  increasing  embankment  height  both  factors  of 
safety  decrease  with  the  ESA  decreasing  at  a faster  rate 
than  for  the  TSA.  Tavenas  et  al.  (1978)  predicted  a rapid 
decrease  in  the  effective  stress  analysis  factor  of  safety 
after  a portion  of  the  clay  foundation  reached  local  yielding 
near  a factor  of  safety  of  1.3  (Figure  4-3).  However,  the 
embankments  they  analyzed  were  not  loaded  to  failure  and 
the  behavior  of  the  Portsmouth  embankment  (Figure  5~3)  does 
not  confirm  this  rapid  decrease  in  the  ESA  factor  of  safety 
near  1.3.  The  safety  factor  from  the  ESA  decreases  at  a 
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rate  faster  than  the  TSA  but  there  does  not  appear  to  be  any 
sudden  change  in  the  slope  to  indicate  local  yielding. 
D'Appoloia  et  al.  (1972)  showed  that  local  yielding  occurs 
at  a factor  of  safety  of  four  for  normally  consolidated  clay 
and  for  heavily  overconsolidated  clay  local  yielding  first 
occurs  at  a factor  of  safety  of  1.8  or  higher.  The  pore 
pressures  measured  a few  hours  after  the  failure  of  the 
Portsmouth  embankment  (Ladd,  1972)  averaged  twenty- five 
percent  higher  than  the  pore  pressures  measured  a few  hours 
before  failure.  These  rapid  pore  pressure  increases  were 
due  to  the  shear  during  failure  and  did  not  provide  any 
indication  of  impending  failure.  The  Interstate-95  embankment 
(Silva-Tulla  et  al.  19'"’6)  also  showed  increased  pore  pressure 
after  failure  but  no  warning  of  impending  failure. 

Of  the  seven  embankments  only  the  New  Liskeard  varved 
clay  and  the  Portsmouth  were  built  to  failure.  Their  factors 
of  safety  were  respectively:  1.00  and  0.90  for  a TSA  and 
1.15  and  0.96  for  the  ESA. 

5.3  LONG  TERM  STABILITY 

The  long  term  stability  of  an  embankment  is  usually 
determined  by  an  effective  stress  analysis  using  hydrostatic 
or  steady  state  pore  pressures.  The  long  term  stability 
determined  by  a total  stress  analysis  is  the  same  as  a "CU" 
class  problem  with  comple+e  consolidation  and  no  pore 
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pressure  dissipation  during  shear.  Although  the  long  term 
class  is  never  the  critical  condition  for  embankments  on 
soft  clays,  a comparison  of  the  factors  of  safety  computed 
from  the  two  types  of  analyses  may  lend  insight  into  their 
use  and  applications. 

The  embankments  considered  are  both  hypothetical* 
examples  using  the  soil  properties  of  two  highly  tested 
clays:  Boston  Blue  Clay  (BBC)  and  Connecticut  Valley  Varved 
Clay  (CWC).  The  BBC  model  embankment  is  described  in 
Appendix  B along  with  the  techniques  used  to  evaluate  the 
factor  of  safety  for  TSA  and  ESA.  The  CWC  model  embankment 
is  described  by  Ladd  and  Foott  (1977)  and  shown  in  Figure  5“5- 

The  end-of-construction  stability  was  determined  for 
a comparison  with  the  long  term  stability.  The  CWC  model 
embankment  had  a "UU"  factor  of  safety  of  1.27  based  on 
corrected  field  vane  (Bjerrum,  1972)  and  Simplified  Bishop 
method  of  slices.  The  BBC  model  had  a "UU"  factor  of  safety 
of  1.02  using  SHANSEP  strength  with  Simplified  Bishop  method 
and  0.94  using  STAB3D  analysis  (Table  5~3) • The  BBC  model 
embankment  was  constructed  slowly  to  allow  increased  strength 
due  to  consolidation. 

For  the  long  term  class  the  ESA  uses  the  drained 
strength  parameters  for  the  soil,  c and  0\  and  static  pore 
pressures.  One-dimensional  consolidation  was  considered  to 
obtain  the  settlement  beneath  the  embankments  but  the  crest 
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elevations  remained  constant.  Both  embankments  were  analyzed 
using  the  Simplified  Bishop  method  of  slices  and  have  factors 
of  safety  of  2.40  and  2.80  for  the  BBC  and  CWC  respectively. 
The  ESA  procedure  is  straight  forward  and  represents  common 
practice  for  determining  the  long  term  stability. 

To  compute  the  long  term  stability  using  a TSA  the 
embankments  were  divided  into  zones  and  the  average  overburden 
of  each  zone  was  used  to  compute  the  average  effective 
consolidation  stress.  For  example,  the  BBC  model  was  divided 
into  three  zones:  centerline  to  crest,  crest  to  toe  and 
beyond  the  toe.  These  zones  also  represent  the  active, 
direct  and  passive  regions  for  shear  strength.  Using  the 
average  effective  stresses  computed  in  this  was  the  SHANSEP 
technique  was  applied  by  computing  the  strength  for  Plane 
Strain  Compression  (PSC),  Direct-Simple  Shear  (DSS)  and 
Plane  Strain  Extension  (PSE)  test  results  for  the  appropriate 
zone.  The  strengths  used  for  both  embankments  were  the 
strengths  on  the  failure  surface  at  failure  (Tff ) . not  the 
maximum  shear  q^,  and  both  were  corrected  for  strain 
compatibility. 

The  CWC  embankment  was  analyzed  by  the  Morgenstern- 
Price  (I965)  method  of  slices  and  the  resulting  minimum 
factor  of  safety  was  1.50.  The  BBC  embankment  was  analyzed 
by  the  Simplified  Bishop  method  of  slices  and  the  resulting 
factor  of  safety  was  1.40.  The  STAB3D  method  of  slices 
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uses  the  anisotropic  strength  properties  and  inclination 
of  the  failure  surface  at  each  slice  to  determine  the  shear 
strength.  The  long  term  factor  of  safety  for  the  BBC  model 
using  the  STAB3D  method  was  1.14.  The  difference  between 
the  Active,  Direct  and  Passive  (ADP)  method  and  the  STAB3D 
method  is  in  the  strength  beneath  the  slope  where  the 'Direct 
Simple  Shear  strength  is  higher  than  the  STAB3D  strength 
(Figure  B-7) • The  factors  of  safety  are  summarized  in  Table 
5-3. 

The  estimation  of  the  SHANSEP  undrained  strengths  may 
be  conservative  because  they  are  based  on  the  K0  stress 
conditions  during  consolidation  and  not  the  in  situ  stress. 
Ladd  and  Edgers  (1972)  showed  the  effects  of  rotating  the 
principle  stress  during  consolidation  on  CAU  DSS  tests. 

The  differences  in  the  two  'types'  of  analyses  for  the  long 
term  class  of  stability  problem  are  still  dramatic  and 
represent  a significant  difference  between  the  two  approaches 
Since  the  long  term  class  is  not  the  critical  class  for 
stability  the  difference  is  not  of  practical  importance 
unless  a rapid  undrained  loading  condition  can  be  expected. 
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TSA  vs  ESA  FACTOR  OF  SAFETY,  ACTUAL  EMBANKMENT  FAILURES 


-3" 

CN 

• 

o 

o 

o 

• 

rH 

-3- 

• 

rH 

o 

SO 

cs- 

co  csj 

• 

• 

csj  • 

o 

o 

• rH 

rH 

SO 

H O 

CS- 

CSJ 

CO 

co 

O 

O os 

Os 

oo 

o 

o 

• 

• • 

• 

• 

• 

• 

tH 

t-h  o 

o 

o 

rH 

rH 

Os 

o 

csj 

'O 

CO 

CO 

Os 

ITS 

o 

oo 

o 

o 

• 

• 

• 

• 

• 

• 

o 

rH 

rH 

rH 

rH 

rH 

Cd 

oo 

CM 

so 

sO 

I csj 

sO 

rH 

oo 

Os 

SO 

O 

COSO 

O 

• 

• 

• 

• 

• 

• • 

• 

H 

o 

H 

rH 

CSJ 

rH  rH 

rH 

>» 

cd  co 

i— i o- 

o os  so 
rH  so 
T3  Os 
CD  -*H 

> TJ 

k TJ  - 

cd  cd  o 

> J (« 


k 0) 
cd  q>  -P 
D CO  C/5 
4<!  to 
to  cd  <8 

•H  O 

►P  cd  o 
►P  X 


X 

• oo 

Z CM  [s_ 
CV  Os 

-OSrH 

rH 

-p 

3 - CD 
O T3  o 
S TO  >J 
to  cd  o 

u 

o 

Ph  cd  ,o 


cs-  - 
Os  Sh 
rH  CD 
0)  > 

<d  - cd 

> T)  <V 
0)  T)  s 
►P  cd 
X «8 

cd 

>>c«  C 
cd  cd 

rH  -p  6 
cd  -p  CH 
ch  o 3 
cd  o cd 
X X 
V 
-p 

< cd  .a 


oo 

csj 

o 

CO 

rH 

o- 

SO 

Os 

O OS 

rH 

•H  O 

• 

§ Os 

T* 

b£)sO 

O 

P Os 

0) 

O r— c 

rP 

CO 

o 

0)  . Cs- 

ctf 

C rH  Os 

S 

•h  cd  r— c 
cd 

S -P  - 

<D  <=8 

f — 1 

cd  c 

cd  >s 

T>  o 

TJ  U 

C33 

to  s~i 

cd  X3  U 

-p  cd 

H cd  o 

-P  X 

■P  P o 

o 

p 

o 

o 

to 

Qu 

1 

-3- 

so 

CSJ 

O 

O 

00 

CSJ 

SO 

CSJ 

O- 

SO 

so 

Os 

Cs- 

SO 

H 

OS£ 

OS 

rH 

rH 

• 

rH 

•k 

•* 

• cd 

e 

u 

co 

•H  * 

0) 

to  -P 

O i-P 

6 csj 

cd  a> 

•H 

rH  CS- 

S 

JC  =8 

cd  os 

cd 

O 

O,  rH 

* rH 

P 

0)  rH 

- Q) 

<=8  - 

M 3 

O >s 

a)  6 

CD  £h 

oo  cd  r 

•P  ^ 3 

> 1 

I X 

Cd  CD  Sh 

a)  cd 

PH  (H 

hOTd  Pi 

x > 

co  rH  cd 

rH 

°8  - 

p,  O ”-3 

— *H 

3 

cd  o « 

SO  CO 

Cs-S  s 

5h 

Os 

Cs- 

o 

1 

1 

w (dp 

M 

m cd  X 

SO 

IN- 

00 

TABLE  5-1 


TSA  vs  ESA  FACTOR  OF  SAFETY 
VARIATION  WITH  EMBANKMENT  HEIGHT 


EMBANKMENT  & REFERENCE  HEIGHT 

FILL 

TSA 

ESA 

1 

Rang  St-George  7 -4m 

Tavenas  et  al.  1978  8.3m 

1.32 

1.25 

2.08 

1.86 

2 

Pang  du  Fleuve  5*  7m 

Tavenas  et  al.  1978  6.7m 

7.8m 

mm 

IIKfl 

2.03 

1.60 

1.35 

3 

Rang  du  la  Concession  4.8m 

Tavenas  et  al.  1978  6.6m 

8.1m 

1.85 

1.66 

1 .40 

2.45 

2.08 

1.75 

4 

Rang  du  Brulfe  6.3m 

Tavenas  et  al.  1978  8.8m 

2.05 

1.55 

2.3 

1.7 

5 

Connecticut  Valley  20.0ft 

Connell  1972  36.7ft 

2.34 

1.46 

2.79 

1.76 

6 

New  Liskeard  Varved  Clay  18. Oft 
Lacasse  & Ladd  1973  20. Oft 

1.20 

1 .00* 

1.28 

1 .05* 

7 

Portsmouth  N.H.  10.0ft 

Boyce  1978  13.0ft 

15.5ft 

18.0ft 

19.5ft 

20.5ft 

21.5ft 

1.88 

1.42 

1.19 

1 .05 
0.98 
0.94 
0.90* 

2.50 

1.88 

1-53 

1.29 

1.14 

1.03 

0.96* 

* Embankment  failure. 


TABLE  5-2 
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TSA  vs  ESA  FACTOR  OF  SAFETY 
MODEL  EMBANKMENTS 


END-OF-CONSTRUCTION* 


METHOD 

BBC 

CWC 

DSS  Simplified  Bishop 
TSA  ADP  Simplified  Bishop 

STAB3D 

FV  Simplified  Bishop 

1.02 

1.02 

0.94 

1.2? 

LONG  TERM 


METHOD 

BBC 

CWC 

ESA  Simplified  Bishop 

2.40 

2.80** 

ADP  Simplified  Bishop 
TSA  Morgenstern-Price 

STAB3D 

1.40 

1.14 

1.50 

BBC  = Boston  Blue  Clay 

CWC  = Connecticut  Valley  Varved  Clay 

DSS  = Direct-Simple  Shear 

ADP  = Active,  Direct  and  Passive  zones. 

* Does  not  include  height  of  fill  required  to 
compensate  for  consolidation  settlements. 

**  Reasonable  failure  surface.  Minimum  occurred 
in  berm  without  passing  thru  clay. 

TABLE  5-3 
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ACTUAL  EMBANKMENT  FAILURES 


TSA  vs  ESA  FACTOR  OF  SAFETY 


ESA 

Factor  of  Safety 


O Field  Vane 
• Corrected  Field  Vane 
■ SHANSEP 

V Unconfined  Compression 
A Triaxial  "UU" 

▲ Average  of  and  A 


* Values  listed  in  Table  5~1 


FIGURE  5-1 


FACTOR  OF  SAFETY 


0.8 


EMBANKMENT  HEIGHT  (FEET) 

* Values  listed  in  Tables  A-4  and  A-5  of  Appendix  A. 

FIGURE  5-3 


AFTER  CONSOLIDATION 


CONNECTICUT  VALLEY  VARVED  CLAY 
MODEL  EMBANKMENT 


FIGURE  5-5 
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CHAPTER  6 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
6.1  SUMMARY  AND  CONCLUSIONS 

This  thesis  considers  the  stability  analysis  of 
embankments  constructed  on  soft  clays.  The  term  'soft'  clay 
implies  that  the  embankment  loading  exceeds  the  preconsolida- 
tion stress  in  the  soil  and  the  clay  becomes  normally 
consolidated.  'Methods'  of  stability  analysis,  'types'  of 
stability  analysis  and  'classes'  of  stability  problems  are 
defined  and  discussed. 

The  'methods'  of  analysis  are  the  techniques  used  to 
determine  the  limiting  equilibrium  factor  of  safety  (Equation 
2-1)  and  involve  the  selection  of  the  critical  failure 
surface  and  the  assumptions  necessary  to  make  statically 
determinate  the  equations  of  equilibrium  (Table  2-1).  Most 
methods  of  analysis  are  reasonably  accurate,  less  than  ten 
percent  difference  between  them,  and  the  choice  of  methods 
is  up  to  the  engineer  based  on  one's  own  experience,  shape 
of  the  assumed  failure  surface  and  technique  (arthimetic, 
graphical  or  computer)  available  for  solution. 


The  ’type'  of  analysis  distinguishes  between  the  choice 
of  expressions  for  calculating  the  shear  strength  of  the  soil. 
The  total  stress  analysis  (TSA)  considers  the  shear  strength 


on  the  failure  surface  at  failure  to  be  independent  of  the 
changes  in  total  normal  stress.  The  strength  may  vary  with 
stress  history,  inclination  of  the  failure  surface  and  level 
of  strain.  The  effective  stress  analysis  (ESA)  uses  the 
Mohr  Coulomb  failure  criteria  (Equation  3-1)  to  express  the 
strength  in  terms  of  the  effective  strength  parameters  *( c 
and  (?)  , the  total  normal  stress  on  the  failure  surface  (<T) 
and  the  pore  pressure  (u). 

The  accuracy  of  either  'type'  of  analysis  is  directly 
related  to  the  ability  to  measure  the  strength  parameters. 

The  TSA  requires  the  single  parameter  of  undrained  strength. 

In  the  past  the  /=0  analysis  was  used  to  obtain  a single 
value  (qf)  for  the  undrained  strength  for  a given  effective 
consolidation  stress  or  natural  water  content.  This  technique 
is  not  considered  accurate  and  in  this  thesis  the  value  of 
the  shear  stress  on  the  failure  surface  at  failure  (T-ff)  with 
consideration  of  the  anisotropic  properties  of  clays  and 
the  effects  of  strain  levels  at  failure,  i.e.  strain 
compatibility,  are  used.  Common  strength  tests  such  as  the 
field  vane,  unconfined  compression  and  "UU"  triaxial 
compression  tests  frequently  do  not  give  accurate  results 
and  thus  an  empirical  correction  factor  for  the  field  vane 
(Bjerrum,  1972)  and/or  more  rational  lab  testing  techniques, 
i.e.  SHANSEP,  which  rely  on  normalized  soil  behavior,  are 
considered  here. 
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The  ESA  strength  parameters,  c and  are  measured  in 
laboratory  tests  and  are  subject  to  numerous  testing  errors 
e.g.  sample  disturbance  and  anisotropy,  as  well  as  test 
interpretation  difficulties  e.g.  use  of  maximum  shear  stress, 
maximum  obliquity  or  line  of  tangency  to  determine  the 
friction  angle.  The  total  normal  stress  is  usually  determined 
by  the  'method'  of  analysis  through  some  assumptions  to 
make  the  equations  of  equilibrium  statically  determinate. 

The  most  difficult  parameter  to  predict  for  the  ESA  during 
the  design  is  the  pore  pressure  which  may  only  be  crudely 
approximated.  During  construction  the  pore  pressures  can 
be  measured  with  piezometers  but  even  then  the  ESA  does  not 
include  the  increased  pore  pressures  generated  during  shear. 

Three  'classes'  of  stability  problems  are  defined, 
based  upon  the  drainage  conditions  during  loading  and  shear 
to  failure.  These  classes  are  labeled  after  similar  lab 
tests:  "UU"  = unconsolidated-undrained  or  end-of-construction, 
"CD"  = consolidated-drained  or  long  term  and  "CU”  = 
consolidated-undrained  or  intermediate  stability  problems. 

In  theory,  each  class  of  stability  could  be  solved  with 
either  type  of  analysis  i.e.  TSA  or  ESA.  However,  in 
practice  on  type  of  analysis  is  usually  better  suited  for  a 
particular  class  of  problem. 

The  TSA  and  ESA  factors  of  safety  are  compared  for 
three  cases:  actual  embankment  failures,  with  variation  in 


n 

1. 
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embankment  height  and  long  term  stability.  The  end-of- 
construction  class  is  the  most  critical  for  embankments  on 
soft  clays  and  are  usually  treated  by  a TSA.  Only  a limited 
number  of  cases  could  be  found  for  actual  embankment  failures 
in  which  both  types  of  analyses  were  used.  These  results 
are  listed  in  Table  5-1  and  plotted  in  Figure  5-1*  Although 
neither  the  TSA  or  the  ESA  appears  very  accurate  in  predicting 
failure  conditions  considering  all  data  equally,  but  when 
only  the  more  accurate  techniques  for  determining  the 
undrained  strength,  i.e.  corrected  field  vane  and  SHANSEP, 
are  used  for  the  comparison  (Figure  5-2) , the  TSA  is  far 
superior  to  the  ESA  for  the  "UU"  class. 

The  factor  of  safety  decreases  with  increasing  embankment 
height  for  both  types  of  analyses  with  both  reaching  a factor 
of  safety  near  1.0  at  failure.  For  the  cases  reported  in  the 
literature  (Table  5-2)  and  the  Portsmouth,  N.H.  embankment 
analyzed  in  Appendix  A,  the  ESA  factor  of  safety  was  higher 
than  the  TSA  at  all  embankment  heights  (Figure  5-4).  The 
ESA  factor  of  safety  decreases  more  rapidly  near  failure 
than  the  TSA  but  there  was  no  sudden  change  as  predicted  by 
Tavenas  et  al.  (1978)  (Figure  4-3)  nor  was  there  any  rapid 
increase  in  piezometer  levels  to  indicate  impending  failure. 
Several  embankments,  Portsmouth  (Ladd,  1972)  and  1-95  (Silva 
Tulla  et  al.  1976),  have  experienced  rapid  increases  in  pore 
pressures  after  failure  due  to  gross  shear  movements  but  in 
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neither  case  were  there  any  indications  of  imminent  stability 
failure.  The  differences  in  the  factors  of  safety  are  due 
to  the  increased  pore  pressures  during  shear  considered  in 
the  TSA  but  neglected  in  the  ESA.  The  difference  does  not 
imply  error  in  either  type  of  analysis  nor  does  it  imply  any 
difference  in  the  stability  of  the  embankment,  but  only  the 
difference  in  strength  expressions  used  in  the  analyses. 

The  long  term  stability  for  both  types  of  analyses  was 
considered  for  two  hypothetical  embankments  constructed  on 
highly  tested  soils:  Connecticut  Valley  Varved  Clay  (Ladd 
and  Foott,  1977)  and  Boston  Blue  Clay  (Appendix  B) . The  ESA 
for  the  long  term  class  uses  the  drained  strength  parameters, 
c and  and  static  or  steady  state  pore  pressures  and  is 
the  approach  most  used  in  practice.  The  TSA  for  the  long 
term  class  is  the  same  as  a "CU"  class  problem  with  complete 
consolidation.  In  both  cases  the  ESA  factor  of  safety  was 
significantly  higher  than  the  TSA  (Table  5-3)  and  again 
demonstrates  the  effects  of  neglecting  the  pore  pressures 
generated  during  shear  to  failure. 

6.2  RECOMMENDATIONS 

The  following  recommendations  are  made  for  the  stability 
analysis  of  embankments  constructed  on  soft  clay. 

The  end-of-construction  class  of  stability  is  the  most 
critical  for  embankments  on  soft  clay.  Design  for  this  class 


of  stability  is  best  conducted  using  a total  stress  analysis 
with  undrained  strength  determined  by  corrected  (Bjerrum,1972) 
field  vane  or  normalized  soil  properties  such  as  SHANSEP 
(Ladd  and  Foott,  197*0*  Strength  based  on  uncorrected  field 
vane  test,  unconfined  compression  or  "UU"  triaxial  tests  are 
subject  to  large  errors  and  may  give  incorrect  results. 

For  field  monitoring  of  the  construction,  particularly 
for  stage  constructed  embankments,  either  type  of  analysis 
may  be  used.  The  TSA  must  be  based  on  normalized  soil 
properties  or  actual  measurements  by  in  situ  testing  or 
sampling  for  each  stage  of  construction.  The  ESA  may  be 
based  on  measured  piezometer  levels  and  the  factor  of  safety 
at  any  instant  of  time  may  be  determined.  The  ESA  is  easier 
to  apply  in  this  case  but  does  not  include  the  effects  of 
pore  pressures  generated  during  shear  and  will  always  give  a 
higher  value  for  the  factor  of  safety  than  the  TSA. 

The  long  term  stability  is  never  the  critical  class  for 
embankments  on  soft  clay  and  therefore  stability  is  not  a 
problem.  The  easiest  type  of  analysis  to  apply  is  the  ESA 
based  on  static  or  steady  state  pore  pressures,  but  the 
results  form  the  ESA  will  be  significantly  higher  than  the 
factor  of  safety  from  a TSA  which  includes  the  pore  pressures 
generated  during  shear. 
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APPENDIX  A 


PORTSMOUTH,  NEW  HAMPSHIRE  TEST  EMBANKMENT 

This  section  analyzes  the  construction  to  failure  of 
an  embankment  located  on  a soft  sensitive  clay.  The 
embankment  was  the  Experimental  Test  Section  of  Inters'tate 
95  in  Portsmouth,  N.H. , which  was  an  extensively  instrumented 
embankment  loaded  to  failure  by  rapid  construction  The 
field  data  are  used  for  two  comparisons:  the  factor  of 
safety  by  total  stress  and  effective  stress  analyses  at 
failure  and  the  variation  in  the  factor  of  safety  with 
embankment  height. 

The  Portsmouth  Experimental  Test  Section  was  selected 
because  of  the  extensive  instrumentation  which  allowed  both 
a total  stress  and  effective  stress  analysis  to  be  conducted. 
This  embankment  failure  was  reported  by  Ladd  (1972)  and 
additional  unpublished  data  were  available  at  M.I.T.  The 
description  of  soil  conditions,  testing  procedures  and 
field  data  are  presented  in  Ladd  (1972)  and  will  not  be 
repeated  here. 

The  soil  properties,  embankment  geometry  and  location 
of  the  failure  surface  are  shown  in  Figures  A-l  and  A-2. 
Because  the  observed  failure  surface  occurred  along  a 
circular  arc,  the  Simplified  Bishop  method  of  slices  (Bishop, 
1955)  was  determined  to  be  a reasonable  method  for  stability 
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analysis  (Wright,  1969) (Wright  et  al.  1969)*  The  stability 
analysis  was  performed  on  LEASE  II  (Dawson,  1972),  a subsystem 
of  ICES,  Integrated  Civil  Engineering  System,  available  at 
M.I.T. 

The  total  stress  analysis  was  performed  on  this  ' end- 
of-construction'  class  problem  using  the  undrained  strength 
determined  by  the  SHANSEP  (Stress  History  And  Normalized 
Soil  Engineering  Properties)  method  with  CK0U  Direct-Simple 
Shear  test  data  and  the  average  preconsolidation  stress 
determined  from  oedometer  tests.  The  strength  values  are 
shown  in  Figure  A-l  and  listed  in  Table  A-2  under  soil  data. 
The  values  of  the  TSA  factor  of  safety  with  variation  in  the 
embankment  height  and  at  failure  are  shown  in  Figure  5-3  and. 
listed  in  Table  A-4.  Run  number  12  was  considered  to  be  the 
most  accurate  and  was  one  percent  lower  than  run  11.  Values 
from  run  11  were  reduced  one  percent  to  reflect  this  decrease 
when  plotted  in  Figure  5~3*  The  factor  of  safety  for  the 
observed  failure  surface  is  the  same  as  reported  by  Ladd 
(1972)  but  the  critical  Bishop  factor  of  safety  is  lower, 

0.90  compared  to  1.01.  This  lower  factor  of  safety  is  due 
to  the  location  of  a different  relative  minimum  and  not  due 
to  any  changes  in  geometry  or  strength.  This  difference 
points  out  the  difficulty  of  evaluating  stability  when 
several  relative  minimums  may  be  calculated. 

The  effective  stress  anaylsis  was  performed  using 


measured  piezometer  readings  throughout  the  foundation  clay 
(Figure  A-2) . The  soil  profile  was  divided  into  layers  and 
a piezometric  line  was  constructed  for  each  layer  based  on 
the  field  readings.  The  layers  included:  a top  layer  of 
sand  with  static  water  table,  a clay  crust  assigned  the 
same  values  of  undrained  strength  as  the  total  stress' 
analysis  due  to  the  difficulties  in  measuring  pore  pressures 
in  this  overconsolidated  crust,  a sand  silt  layer  with 
measured  pore  pressures,  a clay  layer  with  pore  pressures 
averaged  between  the  soft  clay  and  the  sand  silt  layers 
(Figure  A-l ) . The  intermediate  clay  layer  was  selected 
because  of  the  lateral  spreading  of  the  pore  pressures  in 
the  confined  sand  silt  layer.  The  piezometric  line  data 
for  the  layers  described  above  are  given  in  Table  A-3  with 
variation  in  embankment  height.  The  strength  parameters  for 
the  effective  stress  analyses  are  shown  in  Figure  A-l  and 
listed  in  Table  A-2.  The  friction  angle  for  the  soft  clay 
of  21°  corresponds  to  the  value  at  maximum  shear  for  CUC 
tests  and  gives  a factor  of  safety  of  0.96  for  the  critical 
Bishop  surface.  A friction  angle  of  25°  would  be  a more 
reasonable  average  for  this  soft  clay  but  results  in  a factor 
of  safety  of  1.11  for  the  minimum  value,  an  increase  of 
fifteen  percent.  The  factor  of  safety  for  variation  with 
embankment  height  and  at  failure  are  given  in  Table  A-5* 

To  compare  the  strength  along  the  observed  failure 


surface  using  the  two  types  of  analyses,  i.e.  TSA  and  ESA, 
the  embankment  was  divided  into  slices  as  shown  in  Figure 
A— 3 • These  slices  were  drawn  to  intercept  a single  layer 
of  soil  as  used  in  the  total  stress  analysis.  The  stability 
of  this  observed  failure  surface  was  computed  manually  and 
the  average  normal  stress  was  used  to  compute  the  average 
shear  strength  along  the  arc  of  each  slice.  The  value  of 
the  factor  of  safety  for  the  ESA  by  both  LEASE  II  and 
manually  was  1.16  and  for  the  TSA  was  1.01.  The  average 
shear  strengths  along  the  failure  surface  are  computed  in 
Tables  A-6  and  A-7  and  shewn  graphically  in  Figure  A-3« 

The  ESA  gave  higher  values  of  strength  to  the  soil  beneath 
the  crest  and  slope  and  lower  values  of  strength  for  the 
soil  at  the  toe  than  did  the  TSA. 
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INPUT  GEOMETRY 


PORTSMOUTH  N.H.  EMBANKMENT 


POINT  DATA 

No.  X 

(FEET) 

Y 

LINE 

No . 

DATA 

Point 

Left  Right 

Soil# 

1 

0 

-14 

1 

1 

2 

1 

2 

200 

-14 

2 

3 

4 

2 

3 

0 

-10 

3 

5 

6 

3 

4 

200 

-10 

4 

7 

8 

4 

5 

0 

-5 

5 

9 

10 

5 

6 

200 

-5 

6 

11 

12 

6 

7 

0 

-2.5 

7 

13 

14 

7 

8 

200 

-2.5 

8 

15 

16 

8 

9 

0 

0 

9 

17 

18 

9 

10 

200 

0 

10 

18 

19 

9 

11 

0 

5 

ll 

20 

21 

10 

12 

200 

5 

12 

21 

22 

11 

13 

0 

10 

13 

18 

21 

11 

14 

200 

10 

14 

21 

23 

12 

15 

0 

15 

15 

23 

24 

12 

16 

200 

15 

16 

24 

25  • 

13 

17 

0 

20 

17 

25 

26 

13 

18 

50 

20 

18 

26 

27 

13 

19 

200 

20 

19 

27 

28 

13 

20 

0 

22 

20 

28 

29 

13 

21 

58 

22 

21 

29 

30 

12 

22 

200 

22 

22 

30 

31 

12 

23 

102 

33 

23 

31 

32 

12 

24 

107.5 

35-5 

24 

24 

29 

12 

25 

114 

39.5 

26 

115 

41.5 

27 

145 

41 .5 

• 

28 

146 

39.5 

29 

150 

35-5 

30 

158 

33 

31 

176 

33 

32 

200 

30 

TABLE  A-l 
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TABLE  A -3 


TOTAL  STRESS  ANALYSIS  vs  EMBANKMENT  HEIGHT 
PORTSMOUTH  N.H. 


Crest 


Factor  of  Safety 


Failure  Surface 


41.5  1.062 


EFFECTIVE  STRESS  ANALYSIS  vs  EMBANKMENT  HEIGHT 


PORTSMOUTH  N.H. 


CREST 

El. 

(FEET) 

Factor 

of  Safety 

Failure  Surface 



1 

Actual 

Surface 

Critical 

Bishop 

X 

Y 

R 

41.5 

1.156 

.985 

81.25 

45.00 

50.00 

E9 

.995 

85.00 

45.00 

46.66 

K$1 

.963 

80.00 

50.00 

52.56 

U 

40.5 

1.230 

1.0?6 

80.00 

45.00 

50.00 

14 

1.032 

80.00 

51.25 

53.85 

17 

1 .088 

85.00 

40.00 

43.55 

18 

39-5 

1 .341 

1 .141 

81.25 

45.00 

50.00 

15 

1.173 

85.00 

41.25 

44.79 

13 

38.0 

1.519 

1.291 

81.25 

45.00 

50.00 

15 

1.326 

85.00 

41.25 

44.79 

18 

35-5 

2.013 

1-575 

80.00 

45.00 

50.00 

15 

1.531 

80.00 

51.25 

53-80 

18 

33.0 

2.563 

2.017 

80.00 

45.00 

50.00 

15 

1.882 

77.50 

43.75 

46.46 

18 

30.0 

3.504 

2.505 

65.00 

41.25 

48.36 

15 

41.5 

1.316 

1.112 

85.00 

45.00 

46.68 

17* 

RUN  # 

Surface 
Min  (ft) 

Depth 

Max 

ft  Soft 
Clay 

-5 

-14 

21° 

-5 

-14 

21° 

tyjl 

+10 

-14 

21° 

+10 

-14 

25° 

Kui 

0 

-10 

21° 

TABLE  A- 5 
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MANUAL  EFFECTIVE  STRESS  ANALYSIS,  PORTSMOUTH  N.H.  (I) 


SLICE 

1 

0 

2 

AX 

3 

w 

4 

U 

5 

Wsin© 

6 

c*x 

1 

-56.6 

3.3 

1.49 

0 

-1.24 

3.30 

2 

-51.3 

4.0 

4.04 

0 

-3.15 

1.34 

3 

-4o.5 

11.7 

21.64 

26.85 

-14.05 

0 

4 

-28.3 

9.3 

25.83 

35.95 

-12.25 

0 

5 

-18.6 

9-2 

31.96 

42.05 

-10.19 

0 

6 

-9-6 

10.0 

40.38 

47.15 

-6.73 

0 

7 

0 

10.0 

44.28 

50.80 

0 

0 

8 

4-9.6 

10.0 

46.13 

52.75 

+7.69 

0 

9 

18.6 

9.2 

43.19 

53-04 

13-78 

0 

10 

28.3 

9.3 

44.79 

50.82 

21.23 

0 

11 

40.5 

11.7 

49.ll 

44.42 

31.89 

0 

12 

51.3 

4.0 

13.43 

0 

10.48 

1.34 

13 

56 . 6 

3-3 

9.13 

0 

7.62 

3.30 

14 

70.3 

7.7 

10.82 

0 

10.19 

0 

3 


SLICE 

7 

u*x 

8 

W-UAX 

8tanf?  * 

10 

6+9 

11 

M(e) 

F=1.16 

12 

104-11  , 

1 

0 

0 

0 

3.30 

.550 

6.00 

2 

0 

0 

0 

1.34 

.625 

2.14 

3 

19.60 

2.04 

.78 

.78 

.544 

1.43 

4 

20.86 

4.97 

1.91 

1.91 

.723 

2.64 

5 

24.14 

7.82 

3.00 

3.00 

.842 

3.56 

6 

29.42 

10.96 

4.21 

4.21 

.931 

4.52 

7 

31.70 

12.58 

4.83 

4.83 

1.000 

4.83 

8 

32.92 

13.21 

5.07 

5.07 

1.041 

4.87 

9 

30.45 

12.74 

4.89 

4.89 

1.054 

4.64 

10 

29.49 

15.30 

5-87 

5-87 

1.038 

5-66 

11 

32.43 

16.60 

6.40 

6.40 

• 976 

6.56 

12 

0 

0 

0 

1.34 

.625 

2.14 

13 

0 

0 

0 

3.30 

• 550 

5.99 

14 

0 

10.82 

9.08 

9.08 

1.020 

8,89 

g 64.24 


Factor  of  Safety  = 


* 3 = 21' 


64.24 

55-27 


= 1.16 


TABLE  A-6 
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MANUAL  EFFECTIVE  STRESS  ANALYSIS,  PORTSMOUTH  N.H.  (II) 


w-u*x 

M(e) 

N 

n 

su 

1 

1000 

1000 

2 

335 

335 

3 

2.04 

.544 

3-75 

15.39 

94 

245 

4 

4.9? 

.723 

6.87 

10.56 

250 

310 

5 

7.82 

.842 

9.29 

9.71 

367 

355 

6 

10.96 

.931 

11.77 

10.14 

446 

355 

7 

12.58 

1.000 

12.58 

10.00 

483 

355 

8 

13.21 

1.041 

12.69 

10.14 

480 

355 

9 

12.74 

1.054 

12.09 

9.71 

4?8 

355 

15-30 

1.038 

14.74 

IO.56 

536 

310 

16.68 

• 976 

17.09 

15.39 

426 

245 

12 

335 

335 

13 

1000 

1000 

14 

10.82 

1.020 

10.61 

17.31 

514 

514 

Factor  of  Safety 

1 .16 

1.01 

SOIL  PROPERTIES 
PORTSMOUTH  N.H. 


CROSS  SECTION  OF  EMBANKMENT  FAILURE 


PIEZOMETRIC  LINE 


FAILURE  SURFACE 
85,51 

Radius=6l/ 


STRENGTH  ALONG  FAILURE  SURFACE 


PORTSMOUTH  N.H 


SOFT  CLAY 
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HORIZONTAL  SCALE  (FEET) 

FIGURE  A-3 
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APPENDIX  B 

MODEL  EMBANKMENT  ON  BOSTON  BLUE  CLAY 

In  an  attempt  to  compare  the  stability  of  a fully 
consolidated  embankment  by  both  total  stress  and  effective 
stress  analyses,  a hypothetical  embankment  was  placed' on  a 
thirty  foot  layer  of  soil  having  the  properties  of  Boston 
Blue  Clay  (BBC).  The  stress  history  is  assumed  as  shown  in 
Figure  B-l  and  embankment  geometry  as  shown  in  Figure  B-2. 

All  soil  properties  are  those  developed  for  Boston  Blue 
Clay  as  reported  by  Azzouz  (1977)*  The  normalized  strength 
parameters  of  BBC  are  plotted  in  Figure  B-3  for  three 
testing  conditions:  Plane  Strain  Compression  (PSC) , Direct- 
Simple  Shear  (DSS)  and  Plane  Strain  Extension  (PSE).  These 
tests  represent  three  typical  failure  modes  along  a failure 
surface  beneath  an  embankment.  Note  that  the  strengths 
given  in  Figure  B-3  are  the  shear  strengths  on  the  failure 
surface  at  failure  (Tff)  and  not  the  maximum  shear  (q^). 

The  strengths  are  also  corrected  for  strain  compatibility 
because  the  peak  strength  occurs  at  different  strains  for 
the  three  testing  conditions  listed  above  and  hence  the 
peak  strengths  can  not  occur  simultaneously.  The  average 
of  the  peak  strengths  for  the  tests  is  higher  than  the 
average  strength  at  any  given  level  of  strain.  Therefore 
the  values  of  the  Plane  Strain  Extension  and  Compression 
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tests  are  reduced  to  represent  the  strength  mobilized  at  a 
uniform  strain  along  the  failure  surface.  Figure  B-4  shows 
the  undrained  strength  parameters  of  BBC  used  in  the  STAB3D 
method  of  slices.  These  parameters  are  also  for  the  shear 
on  the  failure  surface  at  failure  (Tff)  and  corrected  for 
strain  compatibility.  Figures  B-3  and  B-4  are  tabulated 
in  Table  B-l . 

The  ' end-of-construction'  class  of  stability  was 
evaluated  to  check  that  stability  would  be  a problem  and 
to  compare  the  results  of  Simplified  Bishop  (1955)  and 
STAB3D  (Azzouz,  1977)  methods  of  slices.  The  Simplified 
Bishop  method  of  analysis  used  two  different  strength 
assumptions.  The  first  analysis  employed  the  Direct-Simple 
Shear  (DSS)  strength  along  the  entire  failure  surface  as 
the  average  value  and  the  second  analysis  used  Plane  Strain 
Compression  (PSC) , Direct-Simple  Shear  (DSS)  and  Plane 
Strain  Extension  (PSE)  in  zones  from  the  center  of  the 
embankment  to  the  crest,  the  crest  to  the  toe  and  beyond 
the  toe  respectively  for  an  Active,  Direct  and  Passive  (ADP) 
analysis  (Figure  B-2) . The  STAB3D  program  considers  the 
inclination  of.  the  failure  surface  at  each  slice  to  compute 
the  shear  strength  based  on  the  anisotropic  soil  properties 
input  into  the  program. 

The  normalized  soil  parameters  and  the  stress  history 
of  the  clay  were  used  to  obtain  the  variation  in  the 
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undrained  strength  with  depth,  SHANSEP  approach,  as 
calculated  in  Table  B-3  and  plotted  in  Figure  B-5A.  The 
Simplified  Bishop  method  of  analysis  with  the  DSS  strength 
for  an  average  strength  along  the  entire  failure  surface 
resulted  in  a factor  of  safety  of  1.02  and  failure  surface 
shown  in  Figure  B-6.  Using  the  same  method  but  with  three 
zones  for  strength  also  resulted  in  a factor  of  safety  of 
1.02.  However,  in  the  process  of  finding  the  minimum  factor 
of  safety  the  failure  surface  has  shifted  away  from  the  PSC 
zone  of  higher  strength  and  the  agreement  in  the  factor  of 
safety  is  fortuitous.  The  shift  away  from  the  zone  of 
higher  strength  is  counteracted  by  the  reduction  in  the 
driving  moment  causing  failure  as  a result  of  the  shift. 

The  STAB3D  method  for  the  same  failure  surfaces  obtained 
from  the  Simplified  Bishop  method  resulted  in  factors  of 
safety  about  seven  percent  lower  for  the  end-of-construction 
class  of  stability.  The  major  difference  in  the  two  methods 
is  the  strength  in  the  zone  beneath  the  slope  where  the  DSS 
strength  is  significantly  higher  than  the  strength  determined 
using  the  anisotrophic  properties  in  the  STAB3D  method 
(Figure  B-6). 

The  hypothetical  embankment  was  constructed  at  a rate 
which  allowed  for  an  increase  in  strength  due  to  consolidation. 
Settlements  computed  from  one-dimensional  analysis  were 
included  to  obtain  the  deformed  geometry  after  consolidation 
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(Table  B-2) . The  embankment  crest  was  assumed  to  remain 
level  by  additional  fill  to  compensate  for  the  settlement. 

The  deformed  geometry  shown  in  Figure  B-2  was  used  for  all 
of  the  long  term  analyses. 

The  stability  after  consolidation  was  computed  using 
three  methods:  TSA  by  Simplified  Bishop  with  three  zones 
for  Active,  Direct  and  Passive  (ADP)  analysis;  TSA  by  the 
STAB3D  method;  and  as  ESA  by  Simplified  Bishop  method  of 
slices.  The  effective  stress  analysis  is  the  'type'  of 
analysis  typically  used  for  "CD"  long  term  stability  problems 
with  drained  strength  parameters,  c and  $ , and  hydrostatic 
pore  pressures.  For  this  study  the  value  of  cohesion  c was 
taken  as  zero  and  the  friction  angle  0*  as  thirty  degrees 
which  represents  the  lower  bound  for  BBC  as  reported  by 
Azzouz  (1977)*  The  ESA  factor  of  safety  computed  for  the 
embankment  was  2.40  and  gave  a very  shallow  failure  surface 
(Figure  B-7) • 

The  total  stress  analysis  for  the  long  term  is  the 
same  as  a "CU"  class  stability  problem  with  one-hundred 
percent  consolidation.  The  consolidation  stresses  were 
estimated  using  the  average  overburden  for  the  three  zones 
shown  in  Figure  B-2.  The  zone  from  the  centerline  to  the 
crest  (crest)  consolidated  under  a load  of  1875  PSF  = 15  Ft 
x 115  PCF,  the  zone  from  crest  to  toe  (slope)  consolidated 
under  a load  of  938  PSF  = 7-5  Ft  x 115  PCF,  and  no  load 
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was  considered  beyond  the  toe.  This  assumption  compares 
very  well  with  elastic  solutions  for  embankments  on  a finite 
layer.  Using  the  consolidation  stresses  computed  in  this 
manner,  the  variation  in  undrained  strength  with  depth  was 
calculated  in  Table  B-^  and  B-5  and  is  shown  in  Figure  B-5B. 
The  undrained  strength  calculated  in  this  way  is  conservative 
because  the  rotation  of  the  principal  stress  is  not 
considered  in  determining  the  strength.  Ladd  and  Edgars 
(1972)  used  CAU  DSS  tests  to  show  that  the  rotation  of  the 
principal  stress  caused  an  increase  in  the  strength  compared 
to  CK0 U DSS.  The  total  stress  analysis  factor  of  safety 
should  also  be  smaller  than  for  the  ESA  due  to  the  pore 
pressures  generated  during  shear  to  failure.  The  factor  of 
safety  from  the  ADP  analysis  was  1.40  and  from  the  STAB3D 
analysis  was  1.14  with  failure  surfaces  shown  in  Figure  B-?. 
The  largest  differences  in  the  shear  strength  along  the 
failure  surface  for  the  two  total  stress  analyses  is  again 
in  the  slope  zone  where  the  DSS  strengths  are  higher  than 
those  calculated  from  STAB3D  anisotropic  parameters. 
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